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EFFECT OF VARIABLE VISCOSITY AND THERMAL CONDUCTIVITY ON HIGH-SPEED SLIP FLOW 

BETWEEN CONCENTRIC CYLINDERS 1 

By T. C. Lin and R. E. Steeet 


SUMMARY 

Schamberg was the first to solve the differential equations of 
slip flow, including the Burnett terms , for concentric circular 
cylinders assuming constant coefficients of viscosity and thermal 
conductivity. The problem is solved for variable coefficients of 
viscosity and thermal conductivity in this paper by applying a 
transformation which leads to an iteration method. Starting 
with the solution for constant coefficients, this method enables 
one to approximate the solution for variable coefficients very 
closely after one or two steps. Satisfactory results are shown 
to follow from Schamberg 1 s solution by using his values of the 
constant coefficients multiplied by a constant factor v, leading 
to what are denoted as the effective coefficients of viscosity and 
thermal conductivity. 

INTRODUCTION 

The fact that a gas is not a continuum but actually a 
collection of molecules in rapid but random motion has 
begun to have more and more importance in the aero- 
dynamics of high-speed flow. This is due to the expectation 
that flow through wind tunnels at low pressure or flight at 
extremely high altitudes will not be amenable to analysis 
using classical fluid dynamics. When the mean free path 
of the molecules l is negligible compared with the macro- 
scopic dimension L, which may be wing chord, tunnel 
diameter, and so forth, the classical picture should hold as 
the molecules are so tightly packed together the gas behaves 
just like a mathematical continuum. The ratio IjL is 
defined as the Knudsen number Kn, which is a measure of 
the degree of gas rarefaction. In terms of the better known 
parameters Reynolds Dumber Be and Mach number M, the 
Knudsen number is proportional to MIBe. Hence, although 
not a new parameter, it is a convenient one to use when the 
degree of rarefaction of the gas is of interest. 

Gas dynamics is the continuous-flow regime or Clausius 
gas regime for which the Navier-Stokes equations together 
with the condition of no slip on the boundaries are valid 
and the Knudsen number is extremely small. If the gas 
becomes more rarefied and the Knudsen number increases, 
the effect of slip along the boundaries becomes noticeable, 
although the Navier-Stokes equations remain valid so long 
as the Mach number romains small. This phenomenon has 
been known for over 75 years and has been the subject of 
an extensive study by physicists. Tsien (ref. 1) has sum- 
marized this work very well. During this same period of 


time the solution of Boltzmann’s integral equation by Enskog 
and Chapman, along lines laid down by Hilbert, has led to 
the distribution function for a nonuniform gas as an expan- 
sion in powers of the Knudsen number. This approach 
yields the equations of flow in successive orders of approxi- 
mation, the first order being the Navier-Stokes equations, 
the second order, the Burnett equations, and so forth. The 
third-order approximation has never been carried out and 
the expected complexity of the result does not seem to make 
the attempt worth while, especially as the restrictions on the 
properties of the gas itself are not strictly valid. Chapman 
and Cowling (ref. 2) have presented this theory in their 
well-known treatise. 

Tsien (ref. 1) presented the Burnett equations of motion 
and pointed out that unless the product Mach number times 
Knudsen number (M Kn ) was significant any problem in 
flow could be theoretically solved using the Navier-Stokes 
equations. The question of the proper boundary conditions 
when the higher order Burnett terms are included was 
raised by Tsien but not answered until 2 years later when 
Schamberg, one of Tsien’s students, showed in his doctor’s 
thesis (ref. 3) that the number of boundary conditions 
required for the Burnett equations is the same as for the 
Navier-Stokes equations. 2 While being the same in num- 
ber, the Schamberg boundary conditions are considerably 
more complex, being also expansions in powers of the 
Knudsen number. The first approximations for the slip 
velocity and temperature jump remain essentially the same 
as those used by the physicists in their treatment of low- 
speed slip flows (ref. 4, ch. 8). The second approximation 
which is required when used in conjunction with the Burnett 
equations for high-speed flows is new and, like the Burnett 
terms, of considerable complexity. 

The Burnett equations and the Schamberg boundary 
conditions apply to the domain of high-speed slip flow that 
the aerodynamicist expects to enter first when he leaves the 
domain of classical gas dynamics. Their great complexity 
discourages expectation of a theoretical solution of any 
practically important problem. Hence, the solution of any 
problem, even trivial so far as flows go, is difficult; but, if 
the problem can be set up experimentally, the attempt 
would seem worth while in order to determine the validity 
of the expansions in powers of Knudsen number and a pos- 
sible delineation of the dividing line between gas dynamics 
and slip flow. 


' Supersedes NAOA. TN 2S95, "Effect of Variable Viscosity and Thermal Conductivity on High-Speed Slip Flow Between Ooncentrio Cylinders” by T. O. Lin and R. E. Street, 1953. 
* The correctness of Schamberg's boundary conditions Is not universally acoepted, bnt they are the only ones proposed so far. 


435 



436 


REPORT 1175 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Experimental results have now been obtained for flow 
past spheres (ref. 5) and transverse flow past circular cylin- 
ders (ref. 6) both of which are almost insurmountable 
theoretical problems. Schamberg (ref. 3) solved the plane 
Couette flow problem and the problem of the simple rotation 
between two coaxial cylinders. In order to linearize his 
equations he assumed that the coefficient of viscosity and 
the coefficient of heat conduction of the gas were absolute 
constants. Lin (ref. 7) removed the restriction of constant 
values of these coefficients and recalculated the plane Couette 
flow for a perfect gas with constant specific heats and con- 
stant Prandtl number but with the coefficients of viscosity 
and heat conduction varying as a constant power of the 
absolute temperature. The present investigation does the 
same for the flow between coaxial cylinders. In contrast 
with the plane Couette flow problem the flow between two 
coaxial cylinders rotating relative to each other seems of 
more than academic interest, since an experimental check 
is quite possible and no doubt will be performed in the near 
future. 

The problem is set up in its general form assuming only 
that the Burnett equations and Schamberg boundary condi- 
tions are valid and the flow is steady and stable. Thus the 
streamlines are circles and only the flow in a single plane 
normal to the cylindric axis need be considered. Whether 
such a flow can be stable at high rotary speeds is outside the 
domain of the method used here. A consideration of the 
stability criteria based upon the Navier-Stokes equations 
with slip at the boundary is to be found in reference 8. 
While special laws of dependence of the gas properties are 
assumed, the method is theoretically possible for other laws 
as well as for variable specific heats and Prandtl number. 

This investigation was carried out at the University of 
Washington under the sponsorship and with the financial 
assistance of the National Advisory Committee for Aero- 
nautics. 


FUNDAMENTAL EQUATIONS AND EXPRESSIONS FOR 
STRESS TENSOR AND HEAT-FLUX VECTOR 

It is convenient to start from the general equations of the 
mean motion of a fluid all of whose physical properties vary; 
in Cartesian tensor notation these equations are (ref. 2, 7, 
or 9): 

The continuity equation: 





bxt 


-=0 


( 1 ) 


The momentum equation: 


Du f 


■ LXU '» TP 

P^-pFt- 


bP< 


bxj 


The energy equation: 


DE bu t bg, 

p Dt +Pii tester 0 


( 3 ) 


where p and E are, respectively, the density and the internal 
energy per unit mass; x f is the Cartesian coordinate in the 
physical space; u t , F t , and are, respectively, the com- 
ponents of the velocity of the fluid mean motion, the external 
force per unit mass, and the heat-flux vector in the redirec- 
tion; P tJ is the component of the pressure tensor; and 


D _ b , b 
Dt bt +Uf bx t 


(4) 


is the comoving time derivative or time derivative following 
the -motion as in hydrodynamics. (See appendix A for 
definitions of all symbols.) The summation convention, 
summing over repeated subscripts, is used. 

These general equations cod be derived directly from Max- 
well’s equation of transfer by making use of the properties 
of the summational invariants for molecular encounters 
without determining the form of the molecular-velocity- 
distribution function. The more convenient Cartesian ten- 
sor notation is used rather than the vector-dyadic notation 
preferred by Chapman and Cowling (ref. 2, pp. 61 and 62). 

In terms of the stress tensor r i} , which is defined by 


P i]~p5ij-\rT {j (6) 


(p being the hydrostatic pressure and S fJ , the unit tensor), 
the momentum equation, equation (2), and the energy 
equation, equation (3), take the following forms, respectively: 


Dui 

p Dt' 




DE, but but bq { _ 


( 6 ) 

(7) 


Upon using the continuity equation, equation (1), and the 
first law of thermodynamics 


dQ=dE+p d 


( 8 ) 


together with the definitions dS=dQ/T and H=E-\~ (p/p), 
Q, S, and 3 being the heat received, the entropy, and the 
enthalpy per unit mass of the gas, respectively, it is easy to 
show that 


DE bu { T 

P-m+P —~P T 


bx t 


DS DE Dp 
Dt p Dt Dt 


( 2 ) 


( 0 ) 
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and energy equation (7) becomes 


or 


T DS bu, bg< 

pT Di +Tii ^ } + ^r° 

DE Dp 

p Dt Dt +Tfi bx t + bx. 


( 10 ) 

( 11 ) 


Adding to equation (11) the product of u t and equation 
(0) and considering that r y is a symmetrical tensor, another 
form of the energy equation is obtained: 

p m( R+ l U{U *)~ pUiFt ~^ + ^ i + ^ t (r '^ )== ° (12) 

Equations (3), (7), (10), and (11) are the most familiar 
forms of tbe energy equation. Neglecting terms containing 
the external force F t and setting H equal to c v T, equations 
(6) and (12) yield the general momentum and energy equa- 
tions given by Tsien (ref. 1) and used by Schamberg (ref. 3). 
The sign of r y , being the same as that used by Tsien and 
Schamberg, is opposite to the usual convention. 

It is only through the expressions for the stress tensor r y 
and the heat-flux vector q { that the above momentum and 
energy equations depend on the form of the molecular- 
velocity-distribution function. Let r r y and jp denote the 
rth-order approximations to the stress tensor r y and the 
heat-flux vector g_ { , respectively, and write 


and 


rTtf='Zj Tij™ 
n-0 

(13) 

r2<=Z) 2l <n) 

(14) 


n— 0 


where and are the nth-order corrections to 0 r y and 
0 g ( , respectively. Then the first-order approximation to the 
molecular-velocity-distribution function, that is, the Max- 
wellian distribution, gives (ref. 2, pp. 112, 122, and 123) 


and 


or y =r y (O) =0 

o2<=2< (0) =0 


(15) 

(16) 


which, togethef with equations (6) and (7), yield the Eulerian 
equation of motion 


Dut v , bp n 
p __ pFt+ _ =0 

and the corresponding energy equation 
DE bu< 

p Pf + 2 ? d 7 r 0 


(17) 


( 18 ) 


The second-order approximation to the molecular-velocity- 
distribution function gives (ref. 2, pp. 112, 122, and 123) 


and 


(1) C) bUj 

ir y -r y « = -2 M ^ 

l2l = 2< a>=_ X g 


(19) 


( 20 ) 


where p. and X are the coefficients of viscosity and thermal 
conductivity, respectively, and is the nondivergent sym- 

metrical tensor associated with the tensor that is, 

OXj 

bu { _ 1 /bu t buA 1 bu t 
bz^^Vba^ bxj 3 ba: t w 

In general, any tensor A tj with a bar over it has the following 
meaning (ref. 1): 


A t j — g (AhFA-h ) — — A tt 5 y 


( 21 ) 


Substituting equations (19) and (20) into equations (6) and 
(7) yields the conventional Navier-Stokes equations 

Dut v bp b / bu?v n 

p ~m~ pFi+ ^r 2 b^A M ' (22) 

and the corresponding energy equation of viscous flow 

<*> 

From the definition of the nondivergent symmetrical tensor, 
equation (21), 

b / bR t \_ 1 b / buA lb/ bu , A l b/ buA 

bxj \ bxj) 2 bxj v bx } ) 3 dzt \ bxj 2 dxj \ bxj 

and 

bu { bu t _ bu { bu t 
bxj bxj bxj bxj 

=1 T ftfeV | , fbus biqVH . 

3 L\ba:i bx^J + \ba* bxz) + \br 3 bxj _y 

2L(a5 + W + (5S + W + ls + wJ i0 

Making use of these relations, it is readily seen that equatons 
(22) and (23) check with the momentum and eneigy equations 
of viscous flow (ref. 10). In the present notation the dissipa- 
tion function is simply 
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but bu ( bu t bUi 
$ = 2m 5^ =2m Hxj bx } 


(24) 


which, is always positive and unaffected in form by the fact 
that the coefficient of viscosity p is a variable. 

The third-order correction to the molecular-velocity-dis- 
tribution function, as given by Burnett (ref. 2 or 11), yields 
the second-order corrections to the stress tensor and the 
heat-flux vector for both spherical and Maxwell molecules. 
These corrections are accurate to terms of order (p/p) * * . In 
Cartesian tensor notation they are (refs. 1, 3, and 7) 3 


(S)_7T *** ^ u ‘ 

Ti] 1 p bx k bxj 


K b (p ~l~dp _ bu k duj 2 bu { bu k , 
2 p l_cte A 3 pbxjj dx< bx k bx k .bxj _P 


K* 


m' 


pT dx ( dxj 


_jr M* Z>? I 

ppT bx t bxj 


Z, M ' 


? bTbT 


-u- p 1 bu { bu k 

— ■ "C TT 


8 pT 2 bx t bx/' p bx t dxj 

Maxwell molecules 


(25) 


a m- d 

Qi ~ dl P Tbx, bx t + ' 


dt P T |_3 bx t V bxj +2 bx< bx 

( a M* <$P , n £¥£) Ilji 

V 3 P pb^p UA P bxj + * pTdx,/ dx, 


( 26 ) 


The K ’ s and 0’s are pure constants. Their correct values 
are given below (refs. 7 and 9): 


Maxwell molecules 

*^KKfc)l 


Z 3 ='2 

K 3 =3 

K t =0 

K T dy. 
*~ 6 Jd? 

K .= 8 


Rigid elastic spherical molecules 


r (27) 


K a = 2X1.014 
ir,==3X0.806 
A 4 =0.681 

-%X0.806-0.990 
n al 

K t =8X 0.928 


(28) 


Rigid elastic spherical molecules 


04=3 


(7 T dp\ -] 


15/7 Tdp\ . I 

\2 pdTj 


6l 4 \2 pdT)^ imb 

45 
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AK 

e*=~~X 1.035 

3 

^ (29) 

0 S = — 3X1.030 
04=3X0.806 

(35, T dp\ 
U + pdTj 


II 

CO 

X 

p 

tO 

h- ‘ 
00 

+ 

•fc 1^ 

X 

o 

bo 

o 

05 

Y 



0.150 J 


(30) 


It is noted that the values of 0 a and 0 6 for Maxwellian mole- 
cules given above are different from those given by Chapman 
(ref. 2, pp. 267 to 270). These corrections together with the 
values of 0 for rigid elastic spherical molecules are due to 
Wang Chang and Ublenbeck (ref. 9). For ordinary gases 
(Tjp) (dp/ dT) has a value lying between 1/2 and 1. It 
follows from the corrected expression for 0 S given above 
that all the coefficients K t and 0 t are less than 117/4 instead 
-of 45/4 as given by Chapman and Cowling (ref. 2, p. 270). 

* The last term within the brackets In equation (25) differs from the one given in the refer- 

ences quoted. The correct form for the Burnett terms has been given in reference 9. This 

error was pointed out to the author* by Prof. 8. A. Schaaf and confirmed by Prof. O. A. 
•Trucsden and Mrs. 0. 8. Wang Ohang. 


With the expressions for the stress tensor and the heat- 
flux vector accurate to the second order, the momentum and 
energy equations, equations (6) and (7), become 


p 


Du t 

Dt 





DE , brij 0 bUibUi d 

p Dt ^ bx { ~‘ >l bxj bXi bx 



(31) 


bqT> 


bx { 


— =0 


(32) 


where r tj m and g_ t m are given by equations (25) and (26), 
respectively. Equations (31) and (32) are the momentum 



HIGH-SPEED SLIP FLOW BETWEEN CONCENTRIC CYLINDERS 


439 


and energy equations for slip flow and they reduce to the 
Navior-Stokes equations and the corresponding energy 
equation upon neglecting terms containing r i} m and q f m . 

From equations (19), (20), (25), and (26) it is seen that 
the ratio of a typical term of r« (i) to it u or q ( m to & has the 
same order of magnitude as either pU/pL or p 2 U/\pTL, 
where U and L are the characteristic velocity and length of 
the flow, respectively. Making use of the relations 


and 


B 7 — 1 

y 

a 2 =yBT 


(33) 

(34) 


where R is the gas constant, Cp, the specific heat at constant 
pressure, 7, the ratio of the specific heats, and a, the adiabatic 
speed of sound, together with the definitions that the Rey- 
nolds number Re=pUL/p, Mach number M=U/a, and the 
Prandtl number Pr=<vt/X, it follows that 


ig,= (r— DiW/Be 
From equation (34) and the perfect-gas law 

p=RpT 

it also follows that 


(35) 

(36) 

(37) 


From the kinetic theory of gases (ref. 4, pp. 50 and 147) 

% 

-b- (38) 


and 


p/p=a?/y=irc /8 

p=0A99pcl 


(39) 


where c is the mean molecular speed and l is the mean free 
path of the gas. From equations (38) and (39) 


M V 

0.998 ap 


= 1.2567 l/ Vap 


(40) 


Hence the Knudsen number is 


Kn=l/L=1.256y l,1 M/Re (41) 

Substituting equation (41) into equations (37) and (35), 


^2=0.7967 1,2 MKn 


\pTL 

308505 — 50 29 


■»-™> (V) 


PrMKn 


(42) 

(43) 


For ordinary gases the Prandtl number is approximately 
7 5 

unity and -^7^; therefore pUjpL and /JUfkpTL have the 

same order of magnitude as M Kn. Since the rarefaction of 
the gas increases with the Knudsen number, it is evident that 
for the high-speed flow of a rarefied gas the second-order 
terms t v w and q t m of the stresses and the heat flux become 
relatively important. For 7=1.400 and Pr=0.750, it fol- 
lows from equations (41), (42), and (43) that Kn=lA8QM/ 
Re, pU/pL=0.940M Kn, and p s H/XpPL=O.20lM Kn. 

According to Burnett’s expression for the molecular- 
velocity-distribution function, t w 0) , the third-order correc- 
tions to the stress tensor, will contain terms of the form 
(p/pyCbuiJdxj) 3 , the ratio of which to ir w has the same order 
of magnitude as 


'-M*Kn 2 


(44) 


Similar terms apply for g^/igi. Hence the slip-flow equa- 
tions, equations (31) and (32), cease to be valid if, for a 
given Mach number, the gas is so rarefied that M?Kn 3 is not 
negligible compared with unity. 

The particular problem to be considered in this investiga- 
tion is the slip flow between concentric cylinders. This 
problem has recently been solved by Schamberg (ref. 3, ch. 
VH) for the case of constant coefficients of viscosity and 
thermal conductivity. The present investigation extends 
Schamberg’s solution to include the effect of variable 
coefficients of viscosity and thermal conductivity. 

Assume that the rarefied gas is confined between two con- 
centric cylinders. The inner cylinder, having radius a and 
the uniform temperature T wa , is rotating at constant angular 
velocity its surface velocity being denoted by Z7=<z«„ a ; 
whereas the outer cylinder, having radius b and the unif orm 
temperature T„ t , is held fixed in space. The flow field is 
conveniently described by the cylindrical polar coordinates 
r, <j>, and z, with the z-axis as the axis of the cylinders. 
Assuming that the flow is two-dimensional and steady and 
that the external force can be neglected, 

n,=0, ^=0, ^=0, F { = 0 (45) 

It follows from the symmetry of the problem that 

U.-0, £=0 (46) 

where u, and u r are the velocity components in the z- and 
r-directions, respectively. 

The appropriate equations of motion are obtained by ex- 
pressing the continuity, momentum and energy equations, 
equations (1\ (30, and (321, in plane polar coordinates. 
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This is easily done by making use of the formulas in general 
orthogonal coordinates (refs. 12 and 13). In view of equa- 
tions (45) and (46) it is found that the equation of continuity 
is automatically satisfied. The momentum equations in the 
tangential and radial directions are, respectively (see 
appendix B), 

<48) 

The energy equation becomes 

l»'+?+^(w-?)=° < 49 > 


spaKvXH < f » r 


( 54 ) 

( 56 ) 


dp_pu^ p- 


dr 


+ 


2 I dpi M 3 , M 2 r 2 / duA s 4 / p d?T „ 1 dT 
ZV P Tr\ p+J L~3 {Wj +3 W +2R d?~ B r dP + 




(56) 


In the above equations t t4 , is the viscous shearing stress, % 
and are the normal stresses in the radial and tangential 
directions, respectively, and q T is the radial component of the 
heat-flux vector. 

The explicit expressions for the required components of 
the viscous stress tensor and the heat-flux vector are ob- 
tained by transforming the general expressions, equations 
(25) and (26), respectively, into plane polar coordinates and 
making the reductions required by equations (45) and (46). 
This is done in appendix B. The results are given as follows : 


< 5 °> 

-f [ (n*-f *) ($)’+( !*-§*) ? %- 

i*l 6l)+(n*+l*) (t) 

■»» 

{(n*+t*) @)'+(l*-5*) ? t*+ 

( I rr 2 tv- \ / W\ 8 1 T7- _1_ dT 1 „ rd /I djA 

Vl2^® 3 Ka )\r) Z^pTdr dr^ Ki \_dr\ P dr) 

(52 > 


*'=~ x dF 


(53) 


The expressions for and given in equations (51) and 
(52) are slightly different from those used by Schamberg 
(ref. 3, pp. 152-153). This is due to the correction of the 
last term in the brackets of equation (25) for the Burnett 
terms in the stress tensor. Substituting equations (50), 
(51), (52), and (53) into equations (47), (48), and (49) and 
using equation (27) for the values of the K’s for Maxwell 
molecules give 


T d 

where /3=— Thus one has three equations, equations 

(54), (55), and (56), which together with the perfect-gas 
law, equation (36), are used to find the four dependent 
variables u$, T, p, and p as functions of r. 


BOUNDARY CONDITIONS 

Tt is seen from the previous section that the introduction 
of the higher order approximations to the stresses and the 
heat flux results in an increase in the order of the momentum 
and the energy equations of the fluid mean motion. For 
instance, the first-order approximations to the stresses yr tj , 
equation (19), lead to the conventional Navier-Stokes 
equations, equations (22), of viscous flow, which are partial 
differential equations of the second order; while the second- 
order approximations to the stresses lead to the third-order 
partial differential equations, equations (31), with the 
expression for r,/ 2) given by equation (25) . Since the relative 
importance of the higher order terms of the stresses and the 
heat flux increases with the rarefaction of the gas, this leads 
one to the expectation that the number of boundary condi- 
tions required for the complete evaluation of a slip-flow 
problem should likewise depend on the degree of the rare- 
faction of the gas. However,. it was shown by Schamberg 
(ref. 3) on both physical and mathematical grounds that the 
number of physical boundary conditions required for a slip- 
flow problem is effectively the same as that for the corre- 
sponding flow in the realm of gas dynamics. 

In slip flow, as in gas dynamics, the condition of zero 
relative normal velocity at the boundary still holds, but the 
relative tangential velocity at the boundary is no longer 
zero and the gas temperature differs from the wall tempera- 
ture. These are known as the “slip velocity” and the “tem- 
perature jump,” respectively. 

The expressions for the slip velocity and the temperature 
jump at low Mach number were investigated by Maxwell, 
Millikan, SmoluchowsM, Khudsen, and others (refs. 1 and 
4). If a; and z are the distances tangential and normal to the 
wall, respectively, u and u w are, respectively, the velocity of 
the gas and of the wall in the x-direction, and T and T a are, 
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respectively, the temperature of the gas and of the wall, 
one has from the kinetic theory of gases (ref. 4) the slip 
velocity 


(w)»~ o— ^=0.998 



(57) 


and the temperature jump 

(58) 


where <r is Maxwell’s reflection coefficient, c„ the specific heat 
at constant volume, a, the accommodation coefficient, and l, 
the mean free path. 

For slip flow at high Mach number the above expressions 
give no longer a true description of the physical relations 
and higher order approximations to the slip velocity and the 
temperature jump must be used. A general method for 


the calculation of approximate expressions for the slip 
velocity and the temperature jump, to an arbitrary degree 
of approximation, is given by Schamberg (ref. 3). The 
method applies the laws of conservation of mass, momentum, 
and energy to the infinitesimal layer of gas adjacent to the 
solid surface, referred to as the “sublayer,” and uses the 
nonuniform molecular-velocity-distribution function. 

The first approximations for the slip velocity and the 
temperature jump thus obtained by the use of the first 
approximation to the velocity distribution for a monatomic 
gas of Maxwell molecules agree with the results given by 
previous investigators for low-speed slip flows. 

The second approximations to the slip velocity and the 
temperature jump, which are required in conjunction with 
the second approximations to the viscous stresses and the 
heat flux for high-speed slip flow, are given as follows (ref. 3) : 
The slip velocity 




5 „D fbT\ , R DT bT 9 R Dp bT 4 R DT bp 5 fbiij bw\Z>w~l 
8 Ii Dt\bxJ + 2 T Dt dV 8 p Dt bx/5 p Dt bxj 6 V bz + bxj Dt J 


+ 


(59) 


where j—1 and 2, and the temperature jump 


nww-r.-* if-! §]+(£>' b T (?g+W- i T H i ©+©+ 


where j=l and 2. 


1 /T vmm T dp fbu , , bw\ , nrT1 b*T , n bTbT, „ /d7Y , 1 RT bT bp , RTbpbT 

2 ai( RT ) pbxXbz + b X J +e ^ T bxjbx, +<hB bx j b Xj +esR \bz) + 8 p bxj bxj +e<> p bxj bxj 

1 b 2 T 1 RT bT dpi , /pV ( 1 /DTV , 1 DT D A p\ , 1 1 

14 RT bz* 14 p dz dzj+^p ) j 6 ® T\Dt ) ^HTTt Dt V° g ' T ) +ew(RT > T 1 


■DTbT 
Dt bz 


,-nrr^m D fbT\ bTD n _ , , 1 D 2 T n 1 DTDw 18 bTDw) 

e s(RT) Dt ^ z ) e s {RT) dz Dt [i°g.(pT 1/s )]+2 m {RT) 1 ' 2 Dt Dt 7 bz Dt \ 


(60) 


In equation (60) but not in equation (59) the summation 
convention over two indices is used. This notation has the 
meaning: 

Xl s2x,Xa^y) 

} ( 61 ) 

so that 

bh ( bUj . bw\ bh /bu . bw\ . bh (bv . bw\ 

\bz + b^) (62) 


The time derivative DIDt, when expressed in the above 
notation, is 


D _b , b 
Dt^btr 1 ™ 1 bx } 


(63) 


All of the derivatives in equations (59) and (60) are to be 
evaluated at a point (x, y, z ) as z— >0. These boundary 
conditions are applicable provided that the pressure level 
and the motion of the gas are such that the second approxi- 
mations to the stress tensor 3 t^ and the heat-flux vector 
agi are applicable. This means that the relation M~Kn 2 <^.l 
holds. 
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The values of a u the b’s, c h and the e’s in equations (59) 
and (60) are given as follows (ref. 3) : 4 

6i=— 5.167 


■Uo=% cos <f> 
w a = —u$ sin <j> 
u b = —U 4 , cos <f> 
w t ,=% son 4 > 


} 

I 


(67) 

( 68 ) 


6,=0.8749 



107 

56 


e 7 =— 7.9888 
e 8 = — 5.4912 
e a =— 1.7183 



In general Ci is used instead of Cj. According to the kinetic 
theory for monatomic gases (ref. 2) 7=5/3, X/pc,=5/2, and 
Ci reduces to C\. The values of <r for air vary from 0.79 to 
1.00 while a for air lies between 0.88 and 0.97 (ref. 1). 

To obtain the boundary conditions for the concentric- 
cylinder flow, it is first necessary to express the general 
boundary conditions, as given by equations (59) and (60) 
in the Cartesian coordinate system, in terms of the polar 
coordinate system. Following Schamberg (ref. 3) closely, 
let x and y be the Cartesian axes associated with the polar 
coordinates r and <f> and the auxiliary coordinate systems 
x ai z a and x b , z b by the equations 


x=a+ 2 a =r cos <f> 
y—x a —r sin <f> 
x=b—z b =r cos <f> 
y=—x b =r sin <f> 


} 

} 


(65) 

( 66 ) 


The velocity components u a , w a and u b , w b of the auxiliary 
coordinate systems are related to the tangential velocity u$ as 
shown by figure 1 and equations (67) and (68) : 


* Schamberg’s values for these constants have not been changed to agree with the corrected 
values of the and 0’s given by equations (27) and (29) because the corrections axe negligible 

to the order of approximation used later om 


The partial derivatives with respect to the auxiliary coor- 
dinates are expressed in terms of the partial derivatives with 
respect to r and <f> by means of equations (69) and (70), which 
are easily obtained from equations (65) and (66), respectively: 


£> . ^ b , ,1& 

^--sm^+cos*-^! 

b , b . lb. 
dz a dr r d<j> 


(69) 


b • . & ,1b 

dx b ^ dr r d<j> I 

d d , . 1 d , 

dz b dr r d<f> 


(70) 


All of the first- and second-order partial derivatives appear- 
ing in equations (59) and (60) for the boundary conditions 
can now be transformed into polar coordinates by means of 
equations (67) to (70). After all of the differentiations with 
respect to have been performed, 4> is put equal to zero in 
the resultant expressions, in accordance with figure 1. Be- 
cause of the condition of axial symmetry, equation (46), the 
partial derivatives of u^, T, p, and p with respect to <p all 
vanish. 



The boundary-condition derivatives for both the convex 
(r=o) and the concave (r=b) cylindrical surfaces are given 
in the following table. The transformations for the deriva- 
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tive D/Dt are obtained from equations (63), (45), and (46)-, 
using the fact that the velocity of the wall u a has the value 
U at r=a and zero at r=b, respectively. 



Equivalent radial derivative at — 

Cartesian derivative 

Convex surface 
(r=a) 

Concave surface 
(r=6) 

dT 

dz 

0 

0 

dT 

dz 

dT 

dr 

_dT 

dr 

d*T 

dx* 

1 dT 
r dr 

ldT 
r dr 

&T 

d? 

cPT 

dr 1 

d?T 

dr 1 

&T 
dz dr 

0 

0 

du 

dx 

0 

0 

du 

dz 

du* 

dr 

du$ 

dr 

i Vu 
dz 1 


-m 

d Hi 
dz 1 

dhij, 
dr 1 

<Pu* 

dr» 

d’u 
dz dz 

0 

0 

bio 

dz 

r 

r 

dw 

dz 

0 

0 

d *10 
dz 1 

0 

0 

dPw 

dz 1 

0 

0 

d’to 
dz dz 

d/u^\ 
dr\ r ) 


u 

«♦ 

— v* 

du , dto 
dz""" dz 


4m 

d /du , du>\ 
dz\dz^~ dz/ 

ium 


d /du , duA 
dzVdz” 7 " dz/ 

0 

0 

DT 

Dt 

0 

0 

DPT 

Dl* 

mm 

0 

DfbT\ 
Dt\dx ) 

UdT 
r dr 

0 

Dt\ dz/ 

0 

0 

D /du . du>\ 
Dt\dz + dx) 

0 

0 

Dw 

Dt 

Uu+ 

r 

0 


With these substitutions one obtains from equations (59) 
and (60) the following four boundary conditions: 

(1) At r=a, 

(2) At r= b, 

■(«*)•- 0 -“' ^£[ 40 ?)]+$** <72) 

(3) At r=a, 

Ta=T aa +c 1 (RT a ) u *^ (^) +£z a (73) 

ta \(m / 0 JJ a 

(4) At r= b, 

T,=T„>- Cl <fiT,y'f b (§)+£z, ( 74 ) 

The subscripts a and b denote the evaluation of a particular 
quantity at r=a and r=b, respectively. The quantities 
Xa, X 6 , Z a , and Z 6 are defined by the following equations: 



The values of the numerical constants a u c h e t , e e , and e s are 
given by equations (64). 


SOLUTION OF CONCENTRIC-CYLINDER FLOW 

From the preceding sections it is seen that the problem of 
slip flow between concentric cylinders is reduced to that of 
solving the three differential equations, equations (54), (55), 
and (56), and satisfying the four boundary conditions given 
by equations (71) to (74). For this purpose it is desirable 
to introduce dimensionless constants and variables as follows: 
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H 



(79) 


U 2 

M 2 — — 

“ yRTua 




Xe„=!^ 

fJ-t ca 


*- 

(80) 


p r 




A-wa 

“J 

• 


-4 

* 

li 

© 

p*=p/p a 

x*=x/x tM \ 


u*=uJU 

P*=pjpa 

C P* =C p/ C P m / 

(81; 

T*=T/T m! 

P*=p/Pua 

Rr*=PrjPr w J 



d log, p* 

_ 1 dp* 



dr* 

~P 

* dr* 



■where 

h=b—a (82) 

and c Pica> Hua, and X„ a are the properties of the gas based 
up6n the wall temperature T wa .- Making use of the perfect- 
gas law, equation (36), and of equations (79) to (82), 
equations (54), (55), (56), and (36) become, respectively, 

Ji[fr-)V^.(^)]=0 (83) 


=yM a , 


(u*y- (1 -i)W» 

r*T* ¥{T a *YReu 


< (m*y r 

(^(p^L 


2yM u 


Y^!V 

\dr*J 


+ 


2 M T,iJ ^) +2 


Q Jl fdT*\* 

r* dr* + T*\dr*J' 

T*d log a2 >*1) 
r* dr* Jf 


(l—k) 2 yM w 2 d 
k?p*(T*yRe u 2 dr* 



4 d / r=t d log, p 
3 dr* \ dr* 


d*r* 1 dT* , 2/3 / dT*\* 
(dr*)* r* dr*'T* \dr* ) 

*\ , 2F* dlog,y*l ) 

/H* r* dr* J) 


p*T a *= P *T* (86) 

VELOCITY AND TEMPERATURE DISTRIBUTION 

Integration of equation (83) gives 

u*=Ar*-2 #r*J(^ (87) 

where A and #1? are constants of integration. For p= Con- 
stant, /x* = 1 and equation (87) yields 

u*=jAt*-{~ (87 a) 


. log, r* 


and 


where 


log, 1/& m 


log, r* m=log, 1/k 


r- f 1 * 
Jo M 

i= ra 

^7 Jo M* 


Equations (89) to (91) and (79) insure that 
£=f=0 at r*=l or r =a 


which is Schamberg’s solution. 

From equations (84) and (87) 

(r*)V J* (r*X* gJ)+4(y-l)Pr mi M„ o s B*=0 (88) 

It is convenient to transform the independent variable r* 
into £ or f. The latter are defined, respectively, as 


£=f=l at r*=l/k or r=b 


In terms of £ equation (88) becomes 

u*|(x*^)+4(7- 1) Pr tta M ttO a mW£= 0 
while in terms of the independent variable f it is 


( 86 ) 

(89) 

(90) 

(91) 


■ (92) 
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if 

4 ® _ X r (?- i ) <!f=f - f (93) 

Equation (93) together with equations (90) and (91) gives 
A(0) =A(1) =0. In general A(f)-Cl. 

Integrating once 

~ ^ (F)f+ A <f>df 

and again 

T*=D+mOj ^ Pr^M^nwj £ df X 

f (Vy+W d! (94) 

For constant n and X, ji*=X*=l, ij=l, f=£, and equation 
(94) reduces to 

T*=D+C log. r*~ ( 7 -l)P^M„ a 2 /i 2 (r*)- 2 

which is Schamberg’s solution. 

In terms of f, equation (87) becomes 

(P)f+ A <» dt (95) 

Now let 

(F)A(n =1+ [2mA(f)] n (96) 

n»l 7M 

This series converges for all finite values of mA(f). Making 
use of equation (96), equation (95) gives 

a,*=A+^ [**r+/(T)] (97) 

V 

where 

/(f) =§ ( ~ 2 ^ n+l J ^ [A(f)]"df (98) 

Similarly, from equations (96) and (94), ass uming Pr/c p = 
Constant, 

T*=D+mO>:+^ Pr^M^lPl-J^+giO] (99) 

where 

y(f)=2mj/(f)df (100) 

Now at r=a, f=0, co*=« 0 *, and T*=T*, and, at r=6, 
f = 1, 6 j*=oj 4 *, and T*=T b *. Therefore, from equations (97) 
and (99) 

*>*= Ma *-(co B *-« t *) (!01) 


where 

E(r)=/(f)-i(0) (102) 

and 

T*=T a *-(T a *-T i ,*)f+ 

— 1) ^i^ ( ( t)]r a,6 * )2 [i-f+* 2 f-^+G , (f)] 

(103) 

where 

ff(f)=«G-)-«(o)-ba)-*(o)]r (io4) 

It is noted that F(0) = 0(0) — Q(l) =0. 

Equations (101) and (103) represent the exact solutions 
of differential equations (83) and (84). Since the latter did 
not contain any Burnett terms, it was relatively simple to 
carry out the integration. However, A(f) depends on a 
knowledge of rj and m* which are functions of T* and so 
any numerical solution will reduce to an iteration process 
of approximating A(f) and hence E(f) and (7(f). This is 
also the method of determining the pressure distribution p* 
from equation (85). Before proceeding to p*, it is advan- 
tageous to rewrite the boundary conditions of the section 
"Boundary Conditions” in a similar dimensionless form since 
the constants &>„*, co b *, T a *, and T b * in equations (101) and 
(103) are precisely these boundary values. The actual 
algebraic reduction is carried out in appendix C, leading to 
the following results: 

«.*=l-1.592a 1 Kn'+XJtn* 

(105) 

.y-kXhKn, 2 

(1061 

r.*-i-o.7»e* H ^{r.*-2v*+ 


ai 6 *=1.592ai 


T a *p b * l-P-E(l), 


(7-1) Pr^M u 2 {w a *-u b *y | Knc+ZaKn* 

(107) 

T 6 *=T w6 *+0.796 Cl m \r a *-T b *+ 


(7—1 )Pr ua M w y {u a *-u b y 1 ( 1 ) ^Kn.+ZhKn 2 


Pjltg — y] l\jV Ka 


rs/ 

and X a , X b , Z a , and Z b are the comphcated expressions given 
in equation (C12). Actually, since X a . . . Z b are multi- 
plied by Kn, % above, it is only necessary to use the zero- 
order approximation to these expressions given by equations 
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(C32), (C33), and (C34). The reason for this, of course, is 
that the differential equations and the boundary conditions 
have been derived only to terms in the square of the Khudsen 
number. Thus, it is probably clearer to write equations (105) 
to (108) in the form (eqs. (C27) ) 

6>o * = * (1 + l“a *Kn, + 2 co a *Kn/) 

cc i *= 0 u b *-\- 1 ut,*Kn,+2Ui,*Kn t 1 

T'*=tf*0-+iT,*Kn.+&*Kn.*) 

Tt,*=<jT b *(l^iT b *l£n t -\-iT b *K.n?) 

which explicitly indicates an expansion in powers of Kn, 
The coefficients h bi a *, A co 6 *, h T a *, and h T b * (k= 0, 1, 2) are 
written out in equations (C29) to (C31) of appendix C. 
The subscript written in front of a symbol thus denotes the 
order of the approximation in the expansion. 

PRESSURE DISTRIBUTION 

From equation (85), neglecting the terms containing 

/ m V 

’ which are the order of u> 3 or Kn 3 , the zero-order 
approximation to p* is 


°p*= 


=exp 


(HO) 


The superscript in front of a symbol denotes the order of 
the approximation to the solution of a differential equation. 

From equations (89) and (90) 


Therefore 


r 7j 


'p *=exp M ua 2 («*) 2 dfJ (111) 


In terms of f, equation (115) becomes 
0.6336(1— fc) 2 


or 


>*=exp #] (Hi!) 

Since equation (85) does not explicitly contain terms of order 
Kn, the zero- and the first-order approximation for p* can 
be obtained from equation (110), (111), or (112). From 
equations (112) and (C29) 


.H^).=° 


(113) 


From equations (85), (C5), and (110), log, 3 p*=\og, °p*+ 
a \pKn a 


* K 3 or 


V=(l+o#frOexp (114) 


with exp ( 0 ipKn w f) ^l+o’A Kn a 3 , where 
0.6336(1— &) 2 


o^= 


M fdT*\ 3 } 0.6336(1 — k ) 3 f r ' dr* d 

0 T* 0 \dr*J ) & J i o p*<Zr* 

V op* (3 yMua Xr* dr*) X r*' dr*)' 

r d?T* 1 i fdT^\ 2_i fdry\\ 

oLd(r*) 2 J r* Xdr*r,T* Xdr* ) )) 


du 

Ik* 


(115) 


is obtained from the second-order terms of equation (85) 
upon substitution of the zero-order approximation therein. 


of 


¥ 


fr w? s M - - -f - w*' iwy* •«*** .(f)] «+ ' 

s r ? i ? b= - .(f) .(-f t ^*)+ 2 s w— .(fo- 3 .(?)] « } < n6 > 

Integrating by parts, using equations (C33) and (C34), and neglecting A(f), equation (116) becomes 

"*=-S^iqf (S 5 **) **]f& J] -14 [^]}(w-MWr-8Mi J P>) 4* </f+ 

- 6 /: (*-as*»-*« **) *) 


(117) 
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where 

Ni = (1 -T ut *) (1 - **)*+ ( 7 - 1 )Pr m M m \ 1 
| t VM u *m 

A r a =MJ L m ^(7 — 1 ) Pr„,— | yJ 
Afy=[j y+ (y — l)Pr IOT J M m 3 m 
W*=[| Y+ ( 7 - l)Pr^] 

As the first approximation put 0 m*=o v and assume 


~y 2 =A+A^+A^ 


A=1+(2P-1) 


L-o(p/)- 2 


(1-P) 2 


A oU ^ 0 (P**) 

1_ (1-F) 1 

3 (1-P) S 


Equations (119) and (120) insure that 0 p*=l at f=0 and 
that oP*=oPb* and d u p*/d£=0 at f =1. Then from equation 
(117) 

°^(f) 3 ^= F ° ( f^ 6 )2 [(|^3+AA) (i-^)+ 

(AA+AAKl-^+AAtt-* 8 *)] (121) 


where 


A=(y~1) + (1 - T wb *) ' 


•yk?M a a 


A=[y7-8(7-l)Pr.a]M^ 

A=[f Y-f (y- 1) Pr. a ]M. 
A=[3Y-f (y~ 1) Pruvi] M bs 5 


Then for f = 1 


Y/1\_ 0 - 6336 -M»<lMri0 O/ i \ 7) 1/, ,4N . 

^(l-t-it) 2 1 L 3 7 8 7 

[ i -.w][-? (3+26,; '- i7,; ‘ )+ 

^ Pr ttC ( 3 +2P-47^)+q^Pr <ca ^ — ^ ^ 2 ± F) j | + 

2 jg g ± ^ ( 1 -^)[^] cw> 

The integral in equation (114) remains to be evaluated. 
It is necessary to assume some relation between n* and T* 
which, to avoid too much complexity, will be assumed to be 
the simple relation /x*=P*^ where 0 is a constant whose 
value lies between 0.5 and 1.0 (ref. 4, p. 150). Furthermore, 
in many cases 0«O.9 (ref. 14) so the assumption 0=1 or 
M*= T* will be made in order to simplify the evaluation of 
this integral as well as to enable a more direct expansion of 
all distributions later on. 

Equation (114), with the aid of equations (101) and (96), 
then gives for 0=1 

2 P*=[l+^)Kn 3 \ exp — { ( Ua *- W6 *) s (l-R)- 

4m(o> a *— o> 6 *)[Pa> a *+F(l)co a *— 

[^a> a *+P(l)co 0 *— a> t *] 2 0T 2f — 1) + J(f) } (124) 

where 

</(f) =4m(co a *— w 6 *) s J^ P(f)df— 4m(co n *— &> 6 *)[Fu a *4- 
P(l)co a *-co 6 *]J o r t- 2 fP(f)df+2m(a )fl *- W6 *) 2 X 

r r i-wr)] 2 df+s^^ ft (*.*-«.")*- 

J 0 »=1 nl Jo 

[k*o>*+F(l)<**-<A t ,*]k-t+ (a >a*- W >*)£-1-F(r) Y [A fy)]»df 

(125) 

Upon neglecting all terms containing A(f), the first approxi- 
mation to J (f) is 

<7(D= 1 c7(f)=0 (126) 

while, neglecting terms containing the second and higher 
powers of A(f), the second approximation becomes 

v(r) = (co a *-o,^) 2 j 1 (r)+( ^*~ a,6 * )j 2 (i-)-- 

2(co a *-co s *) (a (127) 


308035— G6 30 
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Ji(r)=^(r)+4mJV(r)dr 
«/,($•) =4ifc*i»*JV* S A (r)df 


(128) 


Equations (124) for the pressure distribution, (101) for 
the velocity distribution, and (103) for the temperature 
distribution in terms of the independent variable f represent 
the solution to the problem, valid to terms in the square of 
the Knudsen number, provided the assumptions of constant 
specific heat, of constant Prandtl number and coefficient of 
viscosity, as well as of heat conductivity, proportional to 
the first power of the absolute temperature, are valid. 

The solution (eq. (124)) for the pressure is more simple 
than Schamberg’s solution in the sense that the present 
assumption p*=T* eliminates the explicit dependence of 
the integral in equation (114) upon T*. Figure 2 shows the 
dependence of 3(1) upon k as given by equation (123) for 
the values of the physical parameters selected for air in the 
section “Case of Air” and discussed there. 



k 


Figure 2. — Burnett terms correction to pressure ratio against 
diameter ratio k for various values of Mach number M m . Pr=0.715; 
7 •= 1.400; 0*=1; a=0.900; <r=l. 


The explicit determination of the expansions of the coeffi- 
cient T) and of the variables £, A(f), F(i), Q( f), and so forth 
in powers of Knudsen number is carried out in appendix D. 

FRICTION COEFFICIENT, SLIP VELOCITY, HEAT TRANSFER, 
TEMPERATURE JUMP, AND PRESSURE RATIO 

FRICTION COEFFICIENT 

It is convenient to define the friction coefficient 0/ as 

G f =-^~ (129) 

The expression for is obtained from equation (50) which, 
upon using equations (81), becomes 

"*—(*#)’ -•*»* » 

From equations (89) and (90) 


* d _y d 

Therefore 


n,U dw* 
Tr ^ am df 

(131) 

where 


v-c-vy-wa 

(132) 

Substituting equation (131) into equation (129) and using 

equations (79), (80), and (81), 


^ % 

(133) 

where 


rjRe,=Re wa 

(134) 


From equations (C21), (C27), and (C29) the value at the wall 
of the inner cylinder becomes 


T>. n 4[1+G“ 0 *— lU^Kn.+ isWa*— iO> b *)Kn/] /10eN 
lie, L fa— £{l+£-[F(l)/(l-£)]} 


Similarly, at the wall of the outer cylinder 


■n n 4&[1+ G^q*~ iUt,*)Kn,+ G^a*~ a u t *)Kn t *l 

e,0fb . l+k-[F(l)/(l-k)] 


(136) 


Since there must be equal torque on both the inner and the 
outer cylinders, (7/»/(7/o=Ar l , as shown by equations (135) 
and (136). If 


then 


C f =,C r (l+ 1 C f Kn.+ 1 C f Kn?) 

4 


Rc t o Gf a 


k{i+k-UFo)i(i-ky\) 


(137) 

(138) 
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Re^Cpr 


4 k 


(139) 


l+i-rofoo/a-*)] 

SUP VELOCITY 

From equations (101), (81), (C28), and (C29), the slip 
velocity at the inner wall is 

u-U 


U 


a *Kn e -\riUa *Kn . 3 = i o> a *Kn , + 2 oi a *Kn 3 (140) 


Similarly, at the outer wall it is 

iU b *Kn t + iu b *Kn/ = (i cc b */k)Kn,+ (nu b */k)Kn. 3 (141) 

where iu a *, iw b *, 3 u> a *, and 2 u 4 * are given by equations (C30) 
and (C31). 

HEAT TKANSFEE 

It is of interest to calculate the heat transfer between the 
rarefied gas and the cylinders. If the dimensionless heat 


transfer g* is defined by 


2 \ T 

tea 


where 


then from equations (53) and (81), 

(1 — k)\* dT* 


2 


ktj 


dr * 


(142) 

(143) 

(144) 


For Prjc p = Constant, \*=n*; so, using equations (89) and 
(90), 

* 1 -k dT* 

2= -5^^r (K5) 

From equation (145) and using equations (C22), (C27), and 
(C29), the dimensionless heat transfer at the inner cylinder is 


l-k?-2m-Q'(<S) 


(iT a *—T wb * 1 T b *)Kn e + ( 3 T a *-T ab * 1 T b *)Kn*+ (y-1) PrM at 

(lWa*— ia t *yKn e a + 2 Gw«*— j 


[1+2 G«a* — i u b *)Kn„-{- 
(146) 


Similarly, that at the outer cylinder is 

2»*=^{ l-T M *+( l T a *-T« b \T 1 ,*)mi.+ (*T*-T ab * i T b *)Kn?+ { y -l)PrM„ 3 (1) [l+2G».*-x»»^n,+ 

G«o*~ ic>3 b *) 3 KnJ l +2 Gu 0 *— ifa b *)Kn?\ j (147) 


If 

g*=oq.*0-+iq.*Kn e + i g*Kn, z ) (148) 

then 

M *=-~~^(l— J 7 6 *+(y— l)Pr M 3 l~^ i ~ 2TO ~~o^ r/ (0) | 

020 km ( 1 * b i)rr ua M aa } 

(149) 

* 1— M, rr * i r i\d u a 1— £*— 2Jehn— 0 O'(l)) 

(150) 

TBMPERATUBE JUMP 

From equation (C27), using equations (C29) and (81), 
the temperature jump at the inner cylinder is 

Ta ^ Tua =iT a * Kn.+aT* Kn 3 (151) 

tea 

Similarly, the temperature jump at the outer cylinder is 

T ^~ T ^ =l T b *Kn,+ t T b * Kn 3 (152) 


The expressions for j T a *, iT b *, /T a *, and 2 T b * are given by 
equations (C30) and (C31). 

PRESSURE RATIO 

The ratio of the hydrostatic pressure at the wall of the 
inner cylinder to that at the outer cylinder is obtained from 
equation (114), upon using equations (121) and (81), as 
follows: 

a-[l+g..‘.f(l)] (163) 

For the case equations (124) and (81) give equation 

(153) the form 

£M1+^.V(1)] exp 2v[1 ^ m j a {[(1-V-Wm)+ 

2{l—k?—2m)F{l) + (k~ 3 - 1)J*(1)J M 3 + 

[-2 (1 -4 s ) +2 (1 +F)m+ (1 -k~ 3 +2m)F(l)] 2« a *« 6 *+ 
(£ _a — P— 4m) (u b *) 3 +J (1) } (154) 

In general tl(l) is a small number and can be replaced by 
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%J{ 1), which is given by equation (D41). From equations 
(154) and (C29) 

o© = 6X5 2^[I=F=^W 

2(1 -P-2 m) oF(l) + {k- 3 -l) o^(l)+gJ(l)] (155) 


For the case of equal wall temperature, T wi *= 1 and, from 
equations (D23) and (C29), 0 F(1)=0. Therefore 

Or in terms of « 


p b \ yM, 
— >=exp 


ua e 


6 0 j? 

(7-DPrAf^ 

120 U 


[( 1 +l<+TO<’+^+|! , + ■■■)- 

0 + S‘' h W‘‘ + • ■ • )] 


(156a) 


From equations (C28) and (154), using equations (81), 
(C29), and (C27) and neglecting the first-order corrections 
to oF(l) and oJ( 1), it follows that 


— 4Pm)+2(l — P — 2m) oF(l)+ 
(A- 2 - 1) oF*(l)](i»«*-§ iij *)+ [-2 (1 -P) +2 (1 +P)m- 
(^_l_2m> 0 F(l)] !«**-!«* !J(1)} (157) 


For the case of equal wall temperature tt F(l)=0 and 

2[— 1+P+ (1 +P)m] m*-± iV *U(1) J (158) 

K 1+£+ ^ ea+ i e3+ f " 4+ • • -) lW6 * + 

(H e+ I^+ • • •)] (168a) 


i Pb 


e iV*(rr—i)P rM a 
80 


lv 


The use of effective values p, and X, and hence of an 
effective Keynolds number Be,, similar to the effective 
Knuds en number Kn, defined by equation (109), is a natural 
consequence of the form of the expressions for the skin- 
friction coefficient, the velocity of slip, the heat transfer, and 
the temperature jump. The use of these effective values 
gives expressions closest in form to those of Schamberg and 
agrees with his when the functions F and O' are zero. Since 
17 as determined in appendix D will, for equal wall tempera- 
tures, increase with the Mach number M m and “curvature” 
1/k as shown in figure 3, so also will n, and X e . 

The limiting values of all expressions when £=1 or e — >0 
agree exactly with the plane Couette flow solution of refer- 
ence 7. 



Figure 3. — Variations of Jij and ov with diameter ratio k for various 
values of Maoh number Af„ a . Pr=0.715; a= 1.400; ^*=1; 7'„!,*=1. 

CASE OF AIR 

Assuming that the gas is air and that the cylinders are 
made of metal such as aluminum or brass, the required 
physical constants are given as follows: 

For the ratio of the specific heats take 

7=1.400 (169) 

a value which is reasonably accurate and has the advantage 
of computational simplicity. It also agrees with the value 
given by the kinetic theory for a diatomic gas whose mole- 
cules have five degrees of freedom. 

The kinetic theory gives, for the Prandtl number, 
Pr =47/(67— 5) =0.737 at 0° C for a diatomic gas -with 
7=1.400 (ref. 4, p. 182). For air, various values have been 
used in the study of the laminar boundary layer in com- 
pressible flow by different investigators. The value Pr= 
0.725 was used by Crocco and Conforto (refs. 15 and 16) 
in 1941; Pr= 0.733, by Brainerd and Emmons (refs. 17 and 
18) in 1942; Pr=0.750, by Schamberg (ref. 3) in 1947; and 
Pr= 0.715, by Cope and Hartree (ref. 14) in 1948. The 

value Pr=l was also used by several •writers (refs. 19 to 21) 
for mathematical simplicity. In the calculations to follow 
the value used will be 

Pr= 0.7 15 (160) 

as suggested by Cope and Hartree in 1948, based on the 
latest data for air given by Kaye and Laby (ref. 22), 1941. 

Maxwell’s reflection coefficient <r is given for air on 
machined brass by K. A. Millikan (see ref. 4, p. 299) as 

cr=1.00 (161) 

The accommodation coefficient a has, according to M. 
Wiedmann (see ref. 1, p. 658), the average value of 

«= 0.900 (162) 0 

for air on metal. This value is also used by Schamberg 
(ref. 3). 

With these physical constants, equations (64) give the 
following constants for slip flow: 
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ai=1.253 ' 
c/= 2.498 
€i=0.251 
. €(=1.911 
« 8 =— 5.491 _ 


(163) 


It is noted that, since air is not a monatomic gas, c/ is used 
instead of ci. 

For the case that ■ 


*_ 1.995* 
lU > “(1 +k)«p b * 


(170) 


The variations of iu a * and x u b * with the ratio * are given 
in figures 5 and 6, respectively. It is noted that j u a * is 
independent of the Mach, number M aa , while x u b * depends 
on Mm, through 0 Pt>*- -As * decreases or the curvature 
increases, it is seen that u a * becomes more negative or the 
slip velocity at the inner cylinder increases in magnitude 
while the slip velocity at the outer cylinder decreases. 

Equations (C29), (166), (159), (160), (163), (D31), and 
(C30) give 


0=1 (164) 

Upon using equations (C29), (159), and (160), equation 
(D16) gives 

(165) 


0.569 f 2 1\ 

10 *(!+*) V l-* 2 m) 


M ua 2 -\- 


0.163 r 1 +F Jr 1 1 M„a* 
*(!+*)(!-**) L4(l-**)m + (l-* s ) 2 2m* J fo 


(171) 


and 


For the case of equal wall temperatures, 


T wb *= 1 
o F(l)=0 


} 


(166) . 


Similarly, using equations (C29), (159), (160), (164) to (166), 
and (D28), equation (D9) gives 


a 0.02045 ri-8 P+* 4 9(1+**) , 31 

oV ~ oV lv (1-F) 2 L 2(l-*v ~4(l-* i )m + ^_ 


(167) 


The variations of both Jx? and o’; for various values of the 
ratio k of the radii of the cylinders at various Mach numbers 
and equal wall temperatures with the ratio of the specific 
heats y= 1.400, Prandtl number Pr=0.715, and the viscosity 
index 0=1 are given in figure 3. 

It is seen that, for *>0.40 and M<4, Ji? and U have 
practically the same value, and therefore equation (165) 
can be used for calculating instead of equation (167). 

At smaller values of * and higher Mach number IvKlv; 
that is, the correction tends to decrease the value of p. 

From equation (156), upon using equations (D41), (159), 
(160), (164), and (166), 


W °- 7Q0Q 

lOg.oPft- qV 


(1— i 4 — 4Pm)— 


ca 

„2 

oV 


^ 02002 ri 2 Pm i-*n - 
(l -* 2 ) 2 1_2 U+ - )+ l-F m J (168) 


The variations of lp b * with the ratio * at various Mach 
numbers and at equal wall temperatures with 7= 1.400, 
Pr=0.716, and 0=1 are given in figure 4. It is seen that 
the pressure ratio \p b * increases rapidly with the curvature 
1 /*, especially at high Mach numbers. 

Upon using equations (C30), (166), and (163), equation 
(140) gives 


\U* 


1.995 

*(!+*) 


Similarly, equation (141) gives 


(169) 


T * 0-569 (l 2 *» \ M aa 2 

1 6 1 +* Vto 1 — **/ o Pb* 

0 A63_ f 1+* 2 ** LI MsL C172I 

(l+*)(l-* 2 )L4(l-* 2 )m + (l-*^) 2 2m 2 J (172) 

The variations of \T* and i T b * with the ratio * for various 
Mach numbers and equal wall temperatures with 7=1.400, 
Pr=0.175, and 0=1 are given in figures 7 and 8, respectively. 
It is seen that as * decreases \T a * becomes much more impor- 
tant than iT b *, the latter going to zero with * for all values 
of M. 



(a) log. Ip b * against h. Dashed curves neglect M* terms. 

Figube 4. — Variations of zero-order pressure ratio with diameter 
ratio * for various values of Mach number. Pr= 0.715; y=l 400- 

0=1; r.»*= 1. 
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Fiqube 6. — First-order slip velocity at inner cylinder against diameter 
ratio k. Pr= 0.715; 7=1.400; 0=1; T„ 4 *=l; <r=l. 


Figure 7. — First-order temperature jump at inner oyllnder against 
diameter ratio k. FV=0.716; 7=1.400; 0=1; <*=0.000. 
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Fiothse 8 . — Firat-order temperature jump at outer cylinder against 
diameter ratio k. Pr=0.716; / y=1.400; |9=1; T,i,*=l; <*=0.900. 



Figure 9. — Firat-order correction to parameter 17 against diameter 
ratio k for various values of Mach number M wa . Pr=0.715; 
•y=1.400; /3= 1 ; T wb *=l; a=0.900; <r=l. 



Figure 10. — Firat-order correction to pressure ratio 1 p b * against di- 
ameter ratio k. Pr= 0.715; T =l-400; 0=1; T wl *= 1; <*=0.900; 
<r=l. ' 


In figure 8 the anomalous behavior of the curves for small 
and large values pf k is due to the large exponential values 
of oPt* for increasing Mach number at all values of k except 
in the neighborhood of k= 1. In fact, for k— 1, 1 T b *= 
0.2845M KO 2 , independent of Q p b *. This is seen upon expand- 
ing equation (172) in powers of 6=1—1?- and letting e— >0. 

From equation (D10), the variations of iv* with the ratio 
k for various Mach numbers and equal wall temperatures 
are found and plotted in figure 9. Since it follows 

from equations (C28) that 17 increases with the Knudsen 
number. 

Making use of equations (D41), (159), (160), and (166), 
equation (158) gives 

„ * n n c\r\ 1 ^ '4Jc*77l iHm i • 

1^6 —0.700 — (i—jrfi ~ (2i w o — iV )~r 


2.800 — 


(1-P) 3 


oV 


o.2°o2 


i r 


(173) 

The variations of i p b * with the ratio k for various Mach 
numbers are given in figure 10. Since ip b *<C.0, it follows 
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that the pressure ratio p b * decreases as the rarefaction of 
the gas, that is, the Knudsen number, increases. 

From equations (135) to (139) and (146) to (150), the 
variations of Be e0 C Ja and Re^G^, x G Sa and i C#, o 2a*, and 
o 2 s* with the ratio k are found and plotted in figures 11, 12, 
13, and 14, respectively. 

The zero-order values of both the skin friction and the 
heat transfer are observed to increase in absolute magnitude 
at the inner and decrease at the outer cylinder with increasing 
curvature. 

Upon using equations (159), (160), and (166), equation 
(123) yields 

o^(l )= MlW Kl-^) + tl- e (p/)“ 2 ][-0.136-1.660P+ 

0.169fc 4 +0.060(l- P)(2+P)/m]} (174) 

The variations of 0 ^(1) with the ratio k for various Mach 
numbers and equal wall temperatures are given in figure 2. 

It is seen that 0 ^(1)]>0, for sufficiently small k, and it in- 
creases rapidly as k decreases. It follows that the effect' of 
the Burnett terms in the differential equation tends to 
increase the pressure ratio p 6 *. 

The coefficient of the last term in equation (123) is plotted 
against k in figure 15. For the case of equal wall temper- 
ature, this term becomes zero. 



Figure 11. — Zero-order friction coefficient at inner and outer cylinders 
against diameter ratio k. Pr-= 0.715; 7=1.400; (3=1; 2 T wt *= 1; 
a=0.900; <7=1. 



Figure 12. — First-order correction to friction coefficient for both inner 
and outer cylinders for various values of diameter ratio k. Pr= 0.716; 
7=1.400; (3=1; T„ t *=l; a=0.900; <r=l. 



Figure 13. — Zero-order heat transfer at inner cylinder against diam- 
eter ratio k. Pr= 0.715; 7=1.400; (3=1; T B i*=l; <*=0.900; 
r=l. 
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Figure 14. — Zero-ordet heat transfer at outer cylinder against diam- 
eter ratio k. Pr<=* 0.716; -y=1.400; >3=1 ; a=0.900‘; 

<r= 1. 


Consider in more detail the special case 


and 

*=0.5 
M ua = 2.0 

From figure 3 

017=1-17 

From figure 4 

Jp 6 *= 2.77 

From figure 5 

l tt a *=— 2.67 

hence 

ico a *=— 2.67 

From figure 6 

1 «„*=0.24 

henco 

i«»*=0.12 

From figure 7 

i!T a *=3.73 

From figure 8 

1 T„*=0.42 

From figure 9 

,ij*=0.96 

From figure 10 

iPt*=— 6.30 

From figure 1 1 

oO/a7?fi« = 5.32 

oOjbRe,—1.2Z 

From figure 12 

iO Ia =- 2.79 

and 

\Cf»= — 2.79 

From figure 13 

020 *= — 1-86 

From figure 14 

o2i,*=0.59 

From figure 2 

o^(l)=4.56 


(175) 

(176) 


(177) 



Figure 15. — Variation of function 0.6336(1— fc)*(2 + i?) /2fc* log, i ivith 
diameter ratio k. 

From equations (C33), (C34), and (113) 



Equation (C32) gives then 


oX d =5A6 

(M?& == — 0.16 

oZa=10.0 

oZ»=-0.14 


( 179 ) 
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a o> a *=17.8 

tu b *=— 0.08 

3 T a *=- 27.0 
,T b *= 0.38 


From equations (D9) to (Dll), using equations (177) and 
(180), 

,7=1.17(1+0^6^.— 5.0ZV) (181) 

Neglecting tlie rarefaction correction to + and F( 1), equa- 
tions (135) and (136) give tlie friction coefficients 

Re e C, a =5.32(l-2.79Kn e +n.9Kn*)) 

> (182) 

Re t O n = 1.33 (1-2.79J&1.+ 17.9Knf) ) 

Similarly, equations (140) and (141) give the slip velocities 


u-U 


= - 2.67JSn»+ 17. S>Kn? 


■+= 0 . 24Kn e —0.16 Kn. 2 


Equations (146) and (147) give the dimensionless heat 
transfer 

q a *= 1.86(1+ 10.5KSi.-83.7jKV) ) 

> (184) 

gj,*= 0.59 (1 — 1.3Kn e -\-8.2Knf) ) 

Equations (151) and (152) give the temperature jumps 


T a -T K 


l ~Z.7ZKn t —27 .OKn? . 


T b — T wb 


Q.42Kn e +0.3 8Kn/ 



Equations (C28) and (154) give the pressure ratio 

2- 6 =2.77(l-6.3Zh,+ . . .) (186) 



Figure 17. — Friction coefficient against effective Knudsen number. 
fV=0.716; 7=1-400; 0=1; T„»*= 1; a=0.900; <r=*l; fc=0.6; 
2 . 



Figure 16. — Parameter 17 against effective Knudsen number. Pr= 
0.715; 7=1.400; 0=1; T wb *=l; a=0.900; <r=l; fc=0.5; Af„=2. 
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Figure 10. — Heat transfer against effective Knudsen number. Pr= 
0.716; 7=1.400; /S=l; T*, i*=l; a=0.900; <r=l; &=0.5; M ma = 2. 



Figure 20. — Temperature jump against effective Knudsen number. 
Pr= 0.716; 7 =1.400; j9=l; T» 6 *= 1; «=0.900; <r=l; fc=0.6; 
M«=2, 


The variations of 17 , Rt,G !a , and Re,C /bl ( u a —U)IU and 
u b /U, q a * and q b *, and {T a —T aa )/T m and (T b —T ab )/T wb with 
the effective Knudsen number Kn, for Pr=0.715, 7=1.400, 
0=1, T ab *= 1, a=0.900, <r=l, *=0.5, and M au =2 are given 


in figures 16, 17, 18, 19, and 20, respectively. Similar re- 
sults could be calculated and plotted for any other choice of 
the parameters, but the enormous amount of work required 
to do so for any range of parameters does not appear to be 
justified at present. This is especially true in that no experi- 
mental data exist for a check. 

The effort of the Burnett (second-order) terms or the Kn, 
terms in equations (181) to (185) is seen to counteract the 
effect of the first-order slip terms or the terms in Kn, in all 
cases except the temperature jump at the outer wall. The 
latter would show the same behavior if l ij and P(l) were not 
neglected. 

SUMMARY OF RESULTS 

A study was made to determine the effects of variable 
viscosity and thermal conductivity on the high-speed slip 
flow between concentric cylinders. The results are sum- 
marized as follows: 

1. Satisfactory estimates of the effect of variable viscosity 
and thermal conductivity upon the velocity and temperature 
distributions were obtained from Schamberg’s solution for 
constant values of these coefficients by basing the friction 
coefficient and the coefficient of heat transfer on the effective 
coefficients fit—vnua and A«=ipw These effective values 
of the coefficients p, and A„ in the case of equal wall tem- 
peratures, increased with the Mach number M and the 
“curvature” 1/k. 

2. Only the expression for the pressure ratio p b /p a was 
significantly different in this case from Schamberg’s solution 
with constant p because of the use of effective values of 
Reynolds number and Knudsen number. 

3. The effect of the Burnett terms in the differential 
equation was more pronounced upon the pressure ratio 
PtlPa- This effect increased with the effective Knudsen 
number Kn„ the Mach number M aa , the temperature 
difference (T Ka —T ub )/T wa , and the curvature 1/k; it was 

measured by the factor o'P of equation (123) as used in 
equation (124) and plotted in figures 2 and 15. The effect 
of the Burnett terms could best be seen by an experimental 
determination of the pressure ratio p b /p a . 

4. The effect of the Burnett terms both in the differential 
equation and in the boundary conditions tended in all cases 
considered to counteract the first-order effect of slip velocity 
and temperature jump on the boundary. 

5. For equal wall temperatures the effect of the Burnett 
terms was to increase the pressure ratio and to increase the 
skin friction on and the heat transfer to both cylinders. 

6. The curvature effect on the behavior of the friction 
coefficient, the slip velocity, and the temperature jump as 
k decreased was such that they all increased at the inner 
cylinder and decreased at the outer cylinder. When k 
approaches unity, the values of all quantities reduce to 
those for plane Couette flow. 


University of Washington, 

Seattle, Wash., December 20, 1951. 
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APPENDIX A 

SYMBOLS 


A 

constant of integration, equation (87) 

Ao, Ai, ... At 

numerical constants, equations (120) 


and (122) 

At ~ Fi p bx { 



tensor 

At) 

nondivergent symmetrical tensor associ- 
* ated with A il} equation (21) 


f adiabatic speed of sound, equation (34) 

a 

< determinant of metric tensor, appendix B 
l radius of inner cylinder 

at 

numerical constant in boundary condi- 


tions, equations (64) 

at) 

fundamental metric tensor, appendix B 

B 

constant of integration, equation (87) 

b 

radius of outer cylinder 

bt, bi 

numerical constants in boundary condi- 


tions, equations (64) 

0 

constant of integration, equation (94) 

'Co, tC\ 

constants of integration, equation (D25) 

C r 

friction coefficient, equation (129) 

c 

mean molecular speed, equation (38) 

Cl, Cl 

numerical constants in boundary condi- 


tions, equations (64) 

Cp 

specific heat at constant pressure 

Cp 

specific heat at constant volume 

V 

constant of integration, equation (94) 

D/Dt 

convective time derivative of hydro- 


dynamics, equation (4) 

V* 

covariant derivative operator, appendix B 

E 

internal energy per unit mass of gas, 


equation (3) 

e U Cj, . . . Cio 

numerical constants in boundary condi- 


tions, equations (64) 

Cy 

symmetric rate of deformation tensor, 


appendix B 

F 

function, equation (102) 

Ft 

component of external force per unit 


mass, equation (2) 

J 

function, equation (98) 

G 

function, equation (104) 

£ 

function, equation (100) 

H 

enthalpy per unit mass of gas, equation 


(9) 

h 

gap between cylinders, b— a 

I 

function, equation (D27) 

J 

function, equation (125) 


K-u K t , . . . Kt 


Kn 

Kn, 

k=ajb 

L 

l 

M 

m=log, j=log e -> 

tC CL 


numerical constants of stress tensor, 
equations (27) and (28) 

Knudsen number, l/L, equation (41) 
effective Knudsen number, equation (100) 

characteristic length of flow 
mean free path of gas, equation (39) 
Mach number, U/a, equation (36) 

equation (89) 


N h N 2> N s , N 4 
P« 

Pr ' 

V 

Q 


it 

Tit 

it™ 

R 

Re 


r 

S 

T 

T a 

t 

U 

U, V, w 


u t 

X, Z 

y, z 

x { 


a 


numerical constants, equations (118) 
component of stress tensor, equation (2) 
Prandtl number, c v p.[\, equation (36) 
static pressure of gas, equation (5) 
heat received per unit mass of gas, 
equation (8) 

component of heat-flux vector, equation 

(3) 

rth-order approximation to 
nth-order correction to <#<, equation (14) 
gas constant 

Keynolds number, equations (36) and 
(80) 

radial distance from center of concentric 
cylinders 

entropy per unit mass of gas, equation 
(9) 

absolute temperature of gas, equation (9) 
absolute temperature of wall 
time, equation (1) 

surface velocity of inner cylinder, au wa 
components of macroscopic velocity in 
x-, y-, and 2-direction, respectively 
component of macroscopic velocity, equa- 
tion (1) 

functions, equations (76) to (78) 
Cartesian coordinates of physical space 
component of Cartesian coordinate of 
physical space 

accommodation coefficient, equation (68) 


T 

A 

St] 


e=l— &*, equation (D18) 


T d 1 

viscosity index, — equation (66) 

Christoffel symbol of second kind, 
pendix B 

ratio of specific heats of gas, c v /c, 
function, equation (93) 
unit tensor, equation (6) 


ap- 
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f 

v 

e 

6i, 0j, . . . 6$ 

X 

M 

i 

p 


Tit 
rT,] 
Ti ,« 
$ 

0^ 

/■V 


variable, equation (90) 
constant, equation (91) 
divergence of velocity, appendix B 
numerical constants in heat-flux vector, 
equations (29) and (30) 
coefficient of thermal conductivity, equa- 
tion (20) 

coefficient of viscosity, equation (19) 
variable, equation (89) 
density of gas, equation (1) 

Maxwell’s reflection coefficient, equation 
(57) 

shearing stress 

component of stress tensor, equation (5) 
rth-order approximation to r tj 
nth-order correction to 0 r y , equation (13) 
dissipation function, equation (24) 
angular coordinate in cylindrical coordi- 
nates 

function, equation (115) 
function, o'PloV 1 , equation (121) 




Subscripts : 

0 , 1 , 2 , ... 

P, i, j, k, l, m 
e 

wa 


Superscripts: 

0 , 1 , 2 , . . . 

“j P, h J, k, l, m 
P 


angular velocity of inner cylinder 
dimensionless angular velocity, equation 
(95) 

numerical constant, equation (C5) 

in front of any symbol, order of approx- 
imation to boundary conditions, equa- 
tions (C27) and (C28) 
covariant vector and tensor indices, 
appendix B 
effective value 

conditions at wall temperature T ua 
indicates covariant differentiation, ap- 
pendix B 

in front of any symbol, order of approx- 
imation to solution of differential 
equation, equation (110) 
contravariant vector and tensor indices, 
appendix B 

exponent in viscosity-temperature 
relation 


APPENDIX B 


TRANSFORMATION OF EQUATIONS OF MOTION TO POLAR COORDINATES 


The equations of fluid mechanics have been written in 
Cartesian tensor form in this report as in references 1 and 3. 
However, the problem considered had to be set up in cylin- 
drical coordinates and Cartesian coordinates were abandoned 
in the process. In order to take into account all possible 
systems of curvilinear coordinates, it is best to express the 
equations of fluid mechanics, including the Burnett terms, 
in general tensor form for any space with a metric form of the 
type (ref. 23) 

<h 2 —a ij dx i dx } 


ordinary partial derivative. The divergence of the velocity 
is the scalar 

u”=j= A (Va «*) 

In other words, take the covariant derivative of the contra- 
variant velocity vector u * and contract the result; a= \a t) \, the 
determinant of the metric tensor a y . By definition the 
covariant derivative is 


By the principle of covariance all that one has to do is to 
express each term in the equations as the proper invariant 
(scalar, vector, or tensor), which reduces to the known form 
when the coordinates are Cartesian. The method and nota- 
tion used are to be found in reference 23. The distinction 
between superscripts and subscripts is necessary. 

There is no difficulty with the equation of continuity. It is 
almost in the form of a scalar equation already. Redefine 
the comoving time derivative 


D_ 

Dt 


dt 


-j-u*D 


* 


where Dt is now the covariant derivative operator. For any 
scalar such as p the co variant derivative is the same as the 



I?-**" 


where the terms rf„ are the Christoffel symbols of the second 
kind. The continuity equation is then 


or 


gf+pe-o 


5f+(»»‘)i,=o 


(Bl) 


for, when covariant differentiation becomes ordinary differ- 
entiation, these expressions reduce to the known equation, 
equation (1) (ref. 23, p. 196). 
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Similarly, the momentum equation is a vector equation 
and the energy equation is a scalar equation. All that is 
needed to give them general tensor form is to place the 
indices properly and to use covariant differentiation. They 
are obviously 


p TF t_pi? * + S s+Tt “ |a=0 (B2) 

P^^+P&+-r afi u a \p+ 2“i«=0 (B3) 

which replace equations (6) and (7) . 

Expansions (13) and (14) remain the same as long as it is 
borne in min d that each ry 00 is a covariant tensor and 
g< to) is a covariant vector. Equations (19) and (20) for the 


second-order approximations remain the same: 

o ^ 

iTy— 2|l-^ 

bT 


‘2<=- X &? 


bUi . 


provided is redefined as the covariant tensor 

^7 = «« =2 ( u * P+%n) — g ©ay 


One must use a tj here as it is the fundamental tensor of the 
space; a {/ is the fundamental contravariant tensor, and 
5/ is the fundamental mixed tensor. 


By continuing this process the expressions for the third-order Burnett terms become 

r {J m =K y ^ Qe tl +K 2 ^ |j^ (Aw+Afli) —^A° la a {J — | +^u\ m u n ^ t a i] — {ei m u m]j +e j m u Mi ) + 1 -f- 

g.t m —d l -^,BT\ { - l r0 2 -^j, (07 , )+2tt t | i T I tJ+O 3 ^7)| t e t 1+^— e* nt+0 s ~T lt A t 


where the covariaht vector Ai—Ft—^ <== ^’ 

and 


&T 


bT 


rp, w yt 

X]ir ~ bx> bx l u bx t 

The Cartesian symbol A i} has been replaced by the covariant 
expression 

Aif=2 (-^y+^iO — 0 A t t (iii 

In order to obtain equations (54) to 56), it is necessary to 
introduce plane polar coordinates r=x 1 and <f>=a ?. The 
metric form is 

ds 2 — dr* -}-r* d <£* 

so a n =<i u =l, Oa i =r 1 =-^> and all other values of Cy are 

(L 

zero while a=r 3 . The only nonvanishing Christoffel symbols 
are 

r|i=-r 

r? 2 =r^=i 

Denote the physical components of the velocity vector by 


u r and u$, where Ur—Ui—u 1 and Similarly, 

the physical components of the vector are q T and q$ and, 
of the stress tensor (ref. 23, pp. 145 to 146), 

Trr=nx = 1 J1 


T,*=- ru=rr ls 
T 


Then 


r 2 2=r 3 T 12 


n 1 i) / ^ . 1 bu$ 
6 

bU, ■ U r ■ 1 bUt 
br r r b<f> 


and the equation of continuity becomes 

If+F ^ (rp ^ )+ r ^ (pw * )=0 


(B4) 


(B6) 


Written out, equations (B2) and (B3) become equations 
(31) and (32) in covariant form, £>*=0.*“ D a , 



Du, 

p ~Ut 
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P~+pe-2^u^ a -D^(XT la )+T^u alfi +^ Ia =0 (B7) 

The components of become, using covariant derivatives 

and the nonvanishing Chris toffel symbols, 

Dui dui, t bu r , biir U 4 bu r V 

'Dt == ~bt +U u '\*=-^+ u '-^+-—-— 


bt 


Du-i bn? , y 

^i*= r 


br 

b 


r b<f> 


buj, 


Dt ~ df" r “ lbi +Uf ^ ^ 


._A£ 


Let Fi=F„ Fi—rFj,, Dip=~> and D 3 p=^- The covariant 
components of e i} are 


en "^ 3 9 


1 / bu T . buj, \ 

SF-^) 


=r (^a 


en=r 


-\-u T 


\b<t> 

while the contravariant ones are 


H 6 ^ 


,u_£r _1 9 

or 3 


s i2 =fi!1= I A A^+I M 

2 Kr 2 b<{> 'r br r 2 / 


Since 


l /but , \ 11 

=?(bf+V— 


3 r 2 


0 


bu r 

%11 -¥ 

bu r 

u w = S*- u * 

bUA 

^ l=r W 

^< 2=r (w +ur ) 

it follows that in equation (B7) the dissipation term is 

>, fbu r \ 2 ,1/1 bu r . but uA 3 , /I but , uA 3 1 2 
+ 2(7 a?+-5F— 7) +h-5*+Fj -3 s 

Consequently, the equations of motion and energy in two- 
dimensional polar coordinates are 


(bUr bUr.Ut bu r ttA w , bp b [~ /du T 1 \1 

p Kw +u ^ + Tb*-— j- pFr+ ^- 2 ^L M fe - 3 e Jr 

1 £> r /I bu r but ■uAl 2/i fbu r 1 but u r \ 
r £><^> l \ r b<j> br r ) J r\br r b<j> r) 


+ 


D a T la ™= 0 

/bu& . but . Ui bUi 

p {-bi + ^ + Tbt- 


m 


U^UA jp ,1 bp 
r ) p ^ + rb^' 


fbE. b'E.Ut bE\, f fbu T 


+ 


AT A bu r , but %\l2_br fl but , u r l.\1 
dr|_ \r b<f> br r) J rd<£[_ \ r b<p r 3 )\ 

T(7^+^-T)+ fl ’ T ■*' , ’- 0 W 

1 fl bu T . but u t\ 2 , A ^ , «rV 1 n2 ~| lb f bT\ 

2 \rb<t> + br r) + \rb<t, + rj 3 G J rdrV X dr/ 

® 10 > 

For 0=0 and p= Constant these reduce to the well-known 
equations for an incompressible fluid (refs. 12 and 13). It 
still remains to find the components of q k m and t*,™ and sub- 
stitute into these equations. For the former the result is 
immediate: 

*- a =°- & e %+* £[f # ( eT >+ 2 w f+i ^ U + 

M 2 f (bUr 1 \ dp 1 bp fl bu r but «#Y| , 

3 P |_Vfr "~3 6 ) ~br + 2? bj> \r ^+~dr~77_| + 

fl fl! i 1 Alr_l but Ur 1 A fill 

r 2 b<f> r 2 3 br 

A f f^h-k 


& P T\_\br 3 
1 bT , 


\ dT, 

__1_ bT { 

CL dz^. . 

but «A1 

/ dP 

2 r d<^> ' 

d</> H 

dr r) J 


7^* 3>[|-r| < 6T >+ 2 0 ^-t) 

2 bT fl but . -mA I | s M 2 [~1/~1 d»r ■ but u Abp 
r b<f> \r b<f> r /J ' 3 pp L2\r dtf> dr r / dr 

1 ^ n f A A^ i M 1 ^1 1 a AQAAAfij. 

3r d</> " l "Vr b<t>'r/r d^J" 1 " " p L2 dr \r b^ 

but uA 1 /I du r | d»j, . 

dr r/^r Vr dflb^ dr r /J - *" 

a A f"i A Aiij. AA l. 

6 pTL2 \r d^ -1- dr r / dr " l_ 

A I ^!_i. ^ 

\r b<f> " r y r d<£ 3r d^> j 


dP 



462 REPORT 1176 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

A longer, more tedious calculation leads to the stress components 


Trr 


18 ) : 


(Du t 

l e 

\ t rr 

p 2 j~2 DAy 


\ dr ~ 

"3® 


V |_3 Dr 

3 r 

1 Du r 
r D4> 


|+|e 

1 

CO| )-< 

)]+x. 


1 ^ j^f2^DT_±DpDT\ +K 

3r 3 D<f> 2 ) +JU f P T \3 dr Dr 3 r 3 d<£ D4>J + s pT |_3 V dr / 3 V 


fDT^_l (\ DT\ t ~ 
r D<f>) _ 


+ 

+ 


—K x p 0 y r ^ + r 3 e J+- K2 p \^ 3 f* 5^3 r 3 dr p |_3 \ Dr r D<t> + r ) \r D<f> r ) 

Du r /Du r 2 \ 4 s (1 Dill . K ^[~1 /l DUf, . Du$ uj\ 2 1 /du r y 2 Du, /I ^ 

^\fr-Z e ) + 9 6 - 2 Vr^ + 7) ] +Kt pll2\r D<j, + Dr r) Q^Dr) ^9 Dr {r D<f> + r J 

2/1 Du* , UrYl.jr m 2 p /2 1 D?T ,21 DT 1 &T\ „ p* /2 1 dp dT 1 dp df\ f jd » T2 /I 
9(r'd^ + 7j dr 3 drV +Xl ppTV3 r 2 d£ d«* 3 dr dr /^'pT^ 1.3 M 

o) m 1 tr M 2 Q fl Dur.Du* «A ~ f fl /l 1 dAA A 2 _2 (1 Du r d^_^\ du r /I 

x^>=r^)= 2 -^i-e^-^+-^- T J+X 2 -|_ 2 ^--^-+--^r; ^ 3 ® \r D<j>‘ Dr r) Dr \r d<y 

, iv ^7j-i^^-r\ + K Jill 

dr \r d<*> + r )y Ks p K r D<j>\Dr r ) +JU p p T 2 \r Dr D<f>^r D<f> Dr J 


dTY_l/dT\n 
D<f> ) 3\. dr ./ J 


\,jr l Rl 
'™ p Tr D(j> dr + 


l^v 0 

6 p 


/I dMr , dw^_w^\ 
\r d7 dr r ) 


Since for the symmetrical case considered in this report Ur = 0 and all partial derivatives with respect to <j> vanish, 
e=0,Ai=— and A=0, it follows that these expressions simplify considerably, becoming simply 

7 n (If 


<Zr®=0 

T «, «« m 2 dT l m* dp fdu* u+\ , 

* <ZA<fr r ) + 

1 p 2 d (du* u<\, P 3 (dUj «A , 

i 0i 7 ^ v^ r ~7; + r 

1 p 2 dT /du* li* 


2 8 pT dr \ dr 


?) 


( 2 )_ 1 jt 1 9 d /i <m 9 /^\ 2 , . o % . 

T " W “8^b? ^~ 2 ^Vp *-J~ 2 UV + U ; +2 r dr J + 

1 k MV 1 ?r ^ r (2 — - — 

Ti^pVdF-T) + 3 E *p B K 2 ~^~r dr)+ 

IrrjLdP dT,2 jP ( dT\* 

3 ^PpT dr 1i + 3 Ks pT K \3F J 


(Bll) 

+ 


1 T r jd_ dpdT_}_ T? tL T?(d 
3 ** p P T dr dr 3 ^'pT \3F. 


(B12) 


r*»=:T* r » = 0 

Equations (B8), (B9), and BIO) reduce for this symmetric 
case to 

-P^+f+^+i (T 7 s, - T ^ 0>) =° 

M&dyA*-?)-* 

But T rr =T rr (2) , T ri , = T Ti , m = — p r^ = r^ <2> , 2r=2r (I) 

= — X^-, and 2« = 2« (2) di this case, so these are precisely 

equations (47) to (49) while equations (B12) and (Bll) 
agree with equations (51) and (52). No second-order terms 
occur in equations (50) and (53) since r r *®==0 and <p <a, =0. 
The second-order term for 2«® does not occur because its 
derivative with respect to <j> is zero in the problem considered. 
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APPENDIX C 

REDUCTION OF BOUNDARY CONDITIONS TO DIMENSIONLESS FORM 


Making use of equations (79) to (81), the boundary con- 
ditions, equations (71) to (74), can be written in dimension- 
less form as follows: 



- 1 

+ ( T a *) 5 


(Cl) 


- , . (j*b*y v 

>" + (p b *y(T a *y Xi " 


„ (T h *y<w r _d_ fu^\i , (p b *y 
1 T tt *p b * L r dr*\r*)i + M\T*y Xr> " 

(C2) 

T*= l+C! m (d^\ “ Za *"* (C3) 


t*-t * f CTv m-, w rn4 x 

1 r a *p 6 * Wr*A + ^ 6 *) 2 (2 7 <1 *) 2Zi ' (G4) 


- R m HwaT a _ (1 —k)y 112 M wa fl— k\ w (n . 

wms* — i — tsz r°- ,9e {—J Zn ™ (C5> 


using equation (41). 

The functions X a *, X t *, Z*, and Z&* are obtained from 
equations (75) to (78) by use of equations (80) and (81) and 
are 
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Upon using equations (89), (90), (91), and (C5), equations 
(Cl) to (C4) become, respectively, 

^ ),g>.+£gv .(07) 


TS- (09) 


T * T *. 

b — -Lwb 

where 

and 


km{T a *) vi Wr/o 
- — 796C ^f-l r/)1 - (^\'Kn.+Z b Kn 2 (OlO) 
Kn,= v Kn ua (Oil) 


« _0.6336(l-A) 2 (^*) s X a * 0.6336(1— &) s f 5 r *dV (. RO , J , 

mn a * /9 t dT*.5dT* 5 , , , ^dto^Y] 

-T (s “ iT +s lif “6 


^ _0.6336(1— i) a 0i 6 *) s X 6 

W(T a *y(p t *y fon\T*Y($ b 

9m^iji*cj* dT* ~| 

8>7 df Jf_i 


* 0.6336 (1 -£) 2 r 5 d* w * ( 6 Ado,* dT* 8 dco* d * . 

^ L“6 T W\ e P) W W Is 1 W d? (Jog ‘ * ; 


^ 0.6336(1 — ky(na*yZa* 0.6336(1— 

£ m 


-ky f u irr ±f du*\ 2 i /dT*Y 1 

vy(T*y ~ Bm\T a *y r l7 tta a vdr/^Vdf/ 14 dr !_t “ 

^ [I 7 m “ t **»* ¥ + (n+‘*>' (t+r "•) ?]},,„ 


^ _0.6336(1-&) 2 Gi s *) 2 Z,* 0.6336(1— 

-Afc = 


tVCr^M 2 


'Pm 2 (T 0 *) 2 (p 


^{^•(f) ! + r,(f)-i«. 


d 2 2’* 
df 2 ' 




[| 


?T*<* 


*%!_1 

dr 


+(^+« 6 ) 


**• ?]L 


( 012 ) 


‘ because of equations (C6), (89), (90), (91), and (95). 

From equations (102) and (98) 

F'(c)=/(r)=f) L ~ 2 ^ +1 ^[A(r)r (oi3) 

F" Qf) -£ (r) — 2mA(f)] (C14) 

n=l fw 

where the prime denotes the derivative with respect to f. 
Equation (93) gives 


A'(f)=£ — 1 
v 


Therefore 


F'(0)=0) 
^(1) = 


,=0 } 

)— o) 

F"(0)=4ra 2 (^-l) 

F"(l)=4jt 2 m 2 ^-l) 


(015) 


(Cl 6) 


(C17) 


and from equations (104) and (100) 

G"(f)=/(f)-ff(l)+ff(0) 

=2m/(f)-p(l)+p(0) (018) 

G"(£)=2mf<£) (019) 

whence 

GT"(0)=<?"(1)=0 (020) 

From equation (101), using equations (Cl6) and (017), 

(yt) 0 = 1 -F-F(l) l_2m) 
(^-)ri-V-F(D ( “ 2Pm) 


/d?oi*\ « a * — «&* 

--A-i-a* 2 -- 


Vdf 2 


TO 


(4m 2 n a * In) 


( df 2 ) 0 l-V-F(l) 
while from equation (103), using equation (020), 


(021) 
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(r— l)Pr wa M ua 2 (u a *— w s *) s r , 7 , „ 

[1 — ^ 2m O' (0)] 


(w),— 

(<W\ _ (7-1) 

\ d? ) 0 [1— Ar 3 — ^(3.)]* 

(#T*\ _ {y~l)Pr u M^ a *-^*Y(AFm -) 

\ d ? )i ' [l-P-i^l)] 2 

( 022 ) 

Putting equations (021) and (022) into equations (07) 
to (CIO) yields 


...♦-1 — 1.5 kll J e ^ m] Kn.+gjZ n > (C23) 


1-k 


(C24) 


r.*=i-o {r.«-r..+ 

(7~ 1 )Pr aa M tca \w*- Ub *y jff } } Kn'+ZjLn? 


465 

Equations (C23) to (C26) give in successive steps the 
zero-order approximation; 


0W o *=l 

ow 6 *=0 

oP 0 + =l 

oT b *=T ab * 

the first-order correction terms: 


(C29) 


lW„ 


* 1.582ai(l— £) 

*[1-A S - 0 P’(1)] 

± 1.592a, (i-lc)l<*(T w >*yi* 


(C30) 


(025) 


r. > =r..«+o.796c,fl- m g)ff‘' 1 |r.._ r ,. + 




(C26) 


It is convenient to express co a *, Wft *, Z a * T**, u„*, p„*, 

and i] in ascending powers of the effective Enudsen number 
Kn,. Thus write 


(027) 


ca a * — ooi a *(l + iw a *Kn t + aO) a *Kn s 2 ) 
u b *=aw b *+iw b *Kn e -{- 3 co b *Kn, s 
7 , a*=o?’ a *(l+iT a *En.+ 2 T a *EO 
T b *=oT b *Q.+ 1 T b *Kn e + i T b *Kn . 1 ) . 

Ua* — tfl-ta * ( 1 H" 1 u a *Kn t + 3 U a *Ejl, Z ) 

«» * = (M * + 1«» *Kn t -j- 2 u 6 *Kn, 2 
Vt>* = oP b *(l + ip 4 *Kn e + 2 p j *KnJ l ) 
y=ov(l+iv*Kn.+2V*Kn t *) 

where the subscript 0, 1, or 2 written in front of a symbol 
denotes the order of the approximation to the boundary 
values. 


(C28) 


U-P-oECL)]^* 

,T *— 0.796g,(l-ft) ( 

“ km 1 1_2 “ 6 + 

( Y _l)p r ][f 2 1— 2m— pg^O) ) 

IT )rr wa M tca [1 _ A 2_ o p (1) j 2 J 

i T * — 0-796^(1—11:) f * i 

1 * moP>*(2U*) 1/2 ( 1 lwi + 

u. ar a 1— 4 s — 2& 2 m— 0 (7'(1)) 

(t l)Pr Ka Mua [ 1 _P_ 0 P( 1 )]2 j 

where oZ(l), oG^O), and oG^OQ are obtained from equations 
(102) and (Cl8) using the zero-order approximation, equation 
(C29), to determine 17 and /**, and the second-order correction, 
terms 

A[l— F— oE(l)] ( l£0a _1 “* ~2 lTa ) +oZ ° 
j, — 1.592a 1 (l k)k 2 (T wb *y ,i * 

2 6 [l-P_oE(l)] 0 p 6 * (ia, “ _lW6 _lT “ + 

liT b *- lPb *^+k 0 X b 

(T-l)Pfil/..=i=g^=^W(2 1 „.*-2 1 „,.- |.(C31> 

§,r. , )j+„z« 

2 T b *= 0-796c!(l— /:) * 1 jV 1— 3T A T * 

mopSiT^*) 111 l 1 '”’ 11 a + 2\} 61 wb ) lIb “ 

(l-T w ^ 1?t *+(T-l)PrM^ 1 -=g^^|^X 

(2 i«o* 2 i«,*+i 1 r # *- 1 r 8 *- lP# *)J + A 
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From equations (C12) and (C29), with n*=T* fi , 


*)mm- 

£ (MogiPZ) jl™YJL /dT^\ ,5 ^ s /doj^ ) 

15 oV ^ A oV dr A o£ 4 oV dr A + o TA1 *~ <A dr Af 


f (032) 


Also equations (C21) and (G22) give, upon using equations (C29), 


fdo}* \ 2m 

o V dr A 1 — — oF(l) 

/dco*\ 2 Pm 

oV'dfA - l-P-oF(l) 

/d?a>*\ 4m 2 

oVdr 2 A~[ 1 -^-qF( 1 )]ot, 

/d?w*\ 4 FmW* 

oVdrV~[l-^-aF(l)] 0 J 


.(?),— <1 - r *** ) 


(7— IjPraJkfm 2 

tl-P-oPCl)] 2 

(7— ljPr^M^ 2 

[l-P-oP(l)] 2 


[l-P-2m- 0 G f, (0)] 

[l-P-2Pm-oG r/ (l)] 


A?P*\ (7— DP^MJ 

oV 3FA~ [l-P-oP(l)] 2 


(4m 2 ) 


/PP*\ ( 7-1 )Pr wa M m * 

oVdr 2 A [l-P-oP(l)] 2 


(4pm 2 ) 


(C33) 


(C34) 


APPENDIX D 


DETERMINATION OF APPROXIMATE 
From equations (90) and (91) 


v=J o M*dr 

(Dl) 

€=- P M*dr 

(D2) 

V Jo 


The kinetic theory of gases (ref. 4) gives 



EXPRESSIONS FOR DISTRIBUTIONS 

(93) yields 


„=/>dr 

(D4) 

P 

V Jo 

(DS) 


(DO) 


tiazT? ) 

that is, > (D3) 

Experiment shows that for most gases the value of 0 lies 
between 0.5 and 1, and it varies slightly with temperature. 
Substituting equation (D3)’into equations (Dl), (D2), and 


It is customary to assume 0 constant and for most gases its 
value is quite close to unity. While the calculations can be 
carried out for any constant value of /S, it is much simpler 
to take 0=1 as is done in the text from equation (124) on. 

Substituting equation (103) into equation (D5) and taking 

0 = 1 , 
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> 1 . ( 77 * 1 . i/'T* T ^ tyPrrgMu ig 2 (t»> 8 * 

* V ( “ r 2 ^ 0 6jr_t " [ 1 — A* — -f(l )] 2 x 

[r^O-W-^+ar)]} (D7) 

where J(f)= 

Substituting equation (103) into equation (D4), with 
(9=1, gives 

7)=i (Z , .*+2 , **) + (7-1)1V 1 bi M i * 1 [i ^g p [i (1+P)- 

CDS 

Upon neglecting the rarefaction correction to T’(l) and 1(1) 
and using equations (C27), (C28), and (C29), this yields as 
successive approximations to 17 

a+r-’H- fcgg^' [| 

(D9) 

0 , l v*=l( l T*+ l T^+2( l <,*- lUi *) (1+**) - 

•i^+o/d)] (DIO) 


017 G^a*+2^6*) + [(l"a*— l"i*) 2 + 


2 w— ,*,»»)] ^^^; [|(i+^ a)] (did 

Now A(0)=A(1)=0 and, for 0<f<l, |A(f)|<g:l. As a first 
approximation assume (superscripts denote order of approxi- 
mation to solutions of the differential equations regardless 
of the boundary conditions) 

A(r)= 1 A(f)=0 (D12) 

It follows from equations (98), ( 100 ), ( 102 ), and (104) that 

y (r) (r ) = x f (r) = x <? (r) =o (D13) 

Substituting these into equations (101), (103), (D7), and 
(D 8 ) gives 

l co*=a> a *- a ‘*--p>* (1-Pt) (D14) 


1 r£if rj-t t (rji * r£ *) | ^1 l) ^^yg da * ^&*) ^ 

(1 &*)* 

(1-P)f-Pq (D15) 

d=| (T a *+T^+ {y ~ 1)P ^^l*~ 03b * y 

(Die) 


N=| (r.*+r»*)+^ (Y-i)Pr M.tiflS-vsr (i+| £+ 

60 ^+40 ^ ) ^ Dl7 ^ 

where 

e=l-P (D18) 

and 


y=^jr a *f-i(r a *-r s *)i- 2 + 

(7-l)PrM w W~-«»*) i j" f 1 (1 w j-gjj 


( 1 — Ar *) 2 


(D19) 


For the second approximation, equation (D 6 ) gives 

which, upon using equations (D15) and (Dl 6 ), yields for 

(3=1 


s A(r)=^{ ( 1 -f)+ 


( 7 — l)PrM wt 2 (a> a *— cjj , *) 1 
(1 -**) 


(D 20 ) 

where ’17 is given by equation (D16). It is noted that 
2 A(0)= 2 A(1)=0. 

Neglecting the second and higher powers of A (f), equation 
(98) yields 

’/(r)=4m 2 JV 2 A(r)df 
which, upon using equation (D 20 ), gives 
! /(i-)=^{ (T a *-T>*) [~^+ (l-m)f+mr 2 ]+ 


( 7—1 )Pr aa M mI 1 (o}a*— 0 } t ,*y 
(l-* 3 ) 

Constant 


— mf(l— f)- 


1 +ft*— 

' 2 d-**) 


+ 


(D 21 ) 


It follows from equation ( 102 ) that 

2 F(f)=i{ (T a *—T t *) [| ( 1 — m) (l-^O+d-m)^!- 

(l- Pf)(P -^)j) 


(D 22 ) 


and 


2(1 -A)* 

2 F(1)=^ (T*-T b *) (i+P-i^T) (D23) 

■z’a)- I A(r.*-r.*)(i+I«+i 2 «*+!«■+ . . .) (D24) 
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where e=l— F as before. 

Substituting equation (D21) into equation (100) yields 
for Pr/c p = Constant 

-y {-(JV-r.*) 

(7-l)PrM^ s (« a *-o> 6 *)Tl-^+2Fm , 

1=F L 2(1-F)m " r 

(1— ^j^p^J|+ t Co+ r Ci f (D25) 
From equations (D25) and (104) 

•<?<r)=4(c r.<-ro {[- 5 ^ 1 +i+* i ]r+ 

(7-1) Pr^V-f.*) 1 ^ 1-F ] 1+F^, [ 

(ss+rS?) (i-w-o—M**- 

(D26) 

The integrals 1(f) defined in equation (D7) for £ and I (1) 
in equation (D8) for i? can now be found from equation 


(D26). They are 

m=£*Gm 


(y— l)PrMaa 2 (u*— t ^*) 2 ( r 1 ~ 4 F , 1 1 . , 
Vl-F) (L 4(l-F)‘ r 2mJ s_r 


_3_ 

2m 


f +4S (tf-l) o- 4 ^- 

(H) ^+5^+165^! 

I(1)= J 0 ‘ 2 <?(r)dr 

(Y-l)PrM„W-«^ 8 [ := |(j= S_ ^ 

J} 


(D27) 


+ 


-F) 2 


KD= 


9(1 +F) 

16(1— F)m 4m 

e 3 


(1328) 


-^[(r a *-p 6 *)(i+£+ . . .)+ 

(y-iyPrM^*-^*) 2 ( 0 +^ £+ . . .)] (3329) 


The derivatives (?'(0) and G' (1) needed in the coefficients X T*, 1 T„*, 2 T a *, and 3 r„* in boundary conditions (C27) are 
also found from equation (D26) as follows: 


{(r.*-r,*) ||l +JS- ' 6 \ J J —+(1 -m) F f +2m (l— m) rF r +2w 1 'i' I /c' t )+ 

( T -l)PrM„.» (».♦-»,♦)■ ri+y i = g + l±g„ W -W f M-rtW--g r ,ir«'1V (D30) 

1 — F L 4 2m 1 — F 1 “ J) 

Thus 

2 G' (0)=i |(P„*-P 6 *) |^2+F-m-^^J+(7-l)-PrM Ka 2 K*-^*) 2 [^fz^+(TZp)2“2m]} (D31) 

2 <?'(0)=2^ [-(P a *-2V) (l+ii e+|j £ 2 + • • -)+^ (y-^PrMJ 1 (* a *- W >*)’ (l+l 6+ • • •)] (3332) 

and 

2 <?'(1)=^ {(Ta*-T»*) [n- 2 F+Fm-^^]+( 7 -l)PrM tra 2 (c a *-co s *) 2 [ 4 ^^+^^“^]} (D33 ) * 
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S6r,(1)= 2ls[ (To * _T6 * ) ( 1+ ^ e+ i fiS+ ■ * 0 + ^ (T_1)Pr ' a/ ““ 2(£Oo *~ CJ6+)a ( 1+ l £+ • • *)] 

Using equation (D22) for *F(f) and equation (D20) for 4 A(f), the integrals J x and J 2 in equations (128) become 
Jo 

=-I \ {(.T a *-T b *) [ 2(1 -m)f+i (3-|) (l-Ft)+(3-m) ^+mr^]+^ ~ 1)Pr |~2^g f _ 

KI + T=f) (l-^)+mf^+(2-m) fi*]} 


0>34) 


(D35) 


J,(D— ^ { (r.*-r,*) g+f f-2»- 5 ^]+ 


(Y- l)PrM w 4 ( Wo *-«^ 4 [^ 


1+F . 2Jchn 


-F) 1 (1-F) 4 m. 


0 (D3<8 


J 2 (f)=4Fm 4 JV 4 r 4 A(f)<2r 

(r.*-2V) [-^(l+“) (*-*-l)+ (l+m)flr*- 

m^]+ (T-l)PrM^^- t o 6 *) 2 |-_2^ 

| (lEp+D^-l) + (2+m) 


( 7 -l )PrM^ (o, a *-a, b *y (± €+ ^ **+ . . .)J (D37) 

Ja(1)= ll^[ (T “* _J ' 6 * ) ( 1 ~I £ “i eJ- nO £3+ • • — l)PrM aa 2 (u a * — w;,*)* ~ (l+^ e+ . . .)] 


From equation (127), when use is made of equations (C29), 

U(t) =o«7i(f) +oJi(f) — 2F _2f JF(f) 


and so 


gj-(l) =M) +^(1) -21F(1) 


JJ(i) “ 24[ (1 - i '-** ) ( 1 + t‘+ • • 


(D38) 

(D39) 

(D40) 


[|a+«+®“ '/]} 

(D41) 

«+i5i2‘'+ • • •)] 

(D42) 
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